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Abstract 


A transient technique was validated for making thermal conductivity measurements. 
The technique incorporated a small, effectively spherical, heat source and temperature 
sensing probe. The actual thermal conductivity measurements lasted 30 seconds. After 
approximately 15 minutes of data reduction, a value for thermal conductivity was obtained. 
The probe yielded local thermal conductivity measurements. Spherical sample volumes less 
than 8cm® were required for the materials tested. Thermal conductivity (and moisture) 


distributions can be measured for relatively dry or wetted samples. 


The technique employs an encapsulated bead thermistor. A thermistor, more commonly 
used as a temperature transducer, has the inherent feature of being readily self-heated. A 
computer-based data acquisition and control system regulates the power supplied to the 
thermistor such that its self-heated temperature response approximates a step change. 
Thermal conductivity is deduced from the transient measurement of the power dissipated by 


the probe as a function of time. 


The technique was used to measure the thermal conductivity of fifteen liquids and five 
insulation materials. Two different thermistors types, glass-encapsulated and 
teflon-encapsulated, were evaluated. Capabilities and limitations of each probe type and the 


measurement technique, in general, were observed. 
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Nomenclature 


R2 


area, m? 
probe effective radius, m 


constant used in equation 16, °C 


kn one ? 
(<=) Alfa: dimensionless 


ATR, 
Vis 
constant pressure specific heat, J/kg-°C 


, -C/W 


i- =. dimensionless 

constant used in equation 16, °C/In(Q) 
feedback error, °C 

current, A 

feedback gain, mA/°C 

thermal conductivity, W/m-°C 

heat transfer rate, W 

power generation per unit volume, W/m? 


thermistor resistance, Q 


statistical parameter 


Nomenclature 


fi = spherical coordinate, m 
T = temperature, °C 

T =T-T,°C 

t = time, s 

V = voltage, V 

y = variable of integration 


Greek Symbols 


a = thermal diffusivity, m?/s 

B = transient heat generation constant, W-s‘/?/m? 

I = steady-state heat generation constant, W/m° 

y = ys 

p = mass density, kg/m? 

o = (= - 1) i. dimensionless 
Subscripts 

a = actual 

b = thermistor bead, theoretical model 

Cc = correction 

d = desired 

f = final 

i = initial 

m = medium 

fe) = output 

i = reference 

Ss = surface 


Nomenclature 


Ss steady-state 


t = thermistor, physical probe 


Nomenclature 
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Introduction 


A relatively inexpensive, portable apparatus for in situ thermal conductivity 
measurements of insulating materials would clearly have a widespread appeal to those 
concerned with energy conservation and management. This study contributes to the evolution 
of such an instrument. Thermal conductivity, especially of insulating materials, is rarely 
measured outside a research or quality control laboratory. A practicing engineer generally 
relies on data tabulated in handbooks or the manufacturer’s literature. Because of the large 
number of materials and product variability, thermal conductivity data are well documented 
for only a relatively small sampling of the commercially available insulation products. 
Moreover, these available values are for new, dry materials rather than for actual installed 
materials under long term operating conditions. Thus, the design engineer must first 
categorize a candidate material and then make an allowance, based on experience, for 
degraded thermal performance resulting from material changes and moisture accumulation. 
Consequently, significant uncertainty exists in estimating heat transfer through typical 
insulations using the tabulated thermal conductivity data. An in situ measuring device would 
provide valuable information on the actual performance of insulation under normal and 


unexpected operating conditions. 


Introduction 1 


To this end, a transient measurement technique was investigated to determine the 
thermal conductivity of insulation materials, especially those materials used in solar energy 
collectors and buildings. The technique incorporates an effectively spherical, small-volume 
probe that is interfaced with a computer-based data acquisition and control system. The 
small-volume probe is both an electrical heating element and a temperature sensor. The 
computer-based data acquisition and control system is used to regulate the probe heating 
element output and to monitor the probe temperature sensor. In practice, after the 
small-volume probe is embedded in an insulation material, it is quickly heated in such a way 
that its temperature response approximates a theoretical step change. The thermal 
conductivity of the insulation material is deduced from knowledge of how the power dissipated 
by the probe varies in maintaining the specified temperature rise over time. Transient 
measurements terminating after 30 seconds yield sufficient information for determining 
thermal conductivity. 

The investigation was divided into key segments. A probe possessing the dual qualities 
of a heat source and a temperature sensor while having essentially a spherical geometry was 
needed. Having found one such probe, an experimental apparatus was configured to simulate 
the analytical model on which the measurement technique is based. The experimental work 
focused upon integrating the probe and the digital hardware and writing the needed software. 
The measurement technique was validated using liquids with well-documented thermal 
conductivities. Liquids with thermal conductivities ranging from 0.1 to 0.6 W/m-°C were used 
because of a lack of standard reference materials in the range of insulation materials (0.03 - 
0.05 W/m-°C). The probe-based measurement apparatus was used to make thermal 


conductivity measurements of five insulating materials. 
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Background 


The generic designation “Small-volume probe” is used even though teflon- and 
glass-encapsulated bead thermistors exclusively comprise the experimentally evaluated 
small-volume probes. This preference to discuss the probe in general terms relates to the 
emphasis given throughout the investigation to building and evaluating the potential of a few 
types of sensors. In a previous investigation, teflon-encapsulated bead thermistors were used 
to make local thermal conductivity measurements of a fibrous insulation batt [1]. Thus when 
this study began, the bead thermistor was viewed as the only probe suited to the 
measurement application. During the investigation, improved probe configurations were 
actively sought and several were evaluated. One other probe design, designated the 
“two-sensor” probe, emerged as worthy of a further, in-depth investigation. The features of 
this probe are described in a later section. As a result of time constraints, however, the 
two-sensor probe was not evaluated experimentally. 

An encapsulated bead thermistor is suited for the small-volume heating and temperature 
measurement application. Basically. a thermistor is a resistance-based sensor conventionally 
employed as a temperature transducer. Thermistors, which are ceramic semi-conductors, 
can be self-heated. Of the different types of thermistors, bead thermistors more nearly 


approximate a small sphere. Encapsulation of the thermistor prevents electrical shorting 


Background 3 


when the probe is submerged in an electrically conducting liquid or embedded in a wetted 


insulating material. 


Thermistor Probe History 


The first Known investigations using the thermistor probe for thermal conductivity 
measurements were performed by researchers in the biomedical field. Chato [2] realized that 
a self-heated bead thermistor could be used to simultaneously measure thermal conductivity 
and thermal inertia (kpC,). Because of an erroneous assumption as to the operation of the 
thermistor, however, Chato documented values of thermal conductivity with an estimated 
error of +20 per cent. Later, Bowman and his associates at MIT [3] built on Chato’s 
preliminary work. The main accomplishment of the MIT researchers was the development 
and solution of a more relevant analytical model of the self-heated thermistor probe 
conducting to a surrounding medium [4,5]. Explicit expressions for two thermophysical 
properties, thermal conductivity and diffusivity, were formulated given the results from one 
test. 

Once the method had been validated, a range of applications for the new probe were 
investigated. Among the MIT investigators, Balasubramaniam and Bowman [6] and Valvano 
et al. [7] quantified the thermophysical properties of different biomaterials. Kravets and 
Larkin [8] used a self-heated thermistor to measure the thermal conductivity and thermal 
diffusivity of foods. Woodbury and Thomas [9] incorporated the self-heated thermistor into an 
investigation to determine indirectly the moisture content of a glass fiber insulation matrix by 


measuring the material’s thermal conductivity. 
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Probe and Application Innovations 


The present work along with that of Woodbury and Thomas and Kravets and Larkin have 
built upon the findings made in the biomedical field. For example, the biomedical work used 
relatively small bead thermistors (nominal diameters of 0.75 mm and 1.10 mm). All of the 
non-biomedical investigations have used the largest, commercially available encapsulated 
bead thermistors (nominal diameters of 2.54 mm and 2.80 mm). Maximizing the probe’s 
Surface area enhances the thermal contact for applications with foods and insulating 
materials. When self-heated, however, the smaller probes used in the biomedical studies 
more rapidly approach the desired step temperature response. For biomaterials as well as 
foods, thermal conductivities ranging from water (0.6 W/m-°C) down to castor oil (0.18 
W/m-°C) make up the desired measurement boundaries. Since the present work and the work 
by Woodbury and Thomas sought to evaluate the thermal conductivity of insulation, the lower 
bound is reduced to 0.03 W/m-°C. 

As compared to the other investigations using the thermistor-based measurement 
technique, the present work was unique in selected areas. For example, the temperature 
response of the thermistor was digitally controlled. Analog controllers were used exclusively 
in all prior investigations. The change to an all-digital system permitted a more efficient data 
acquisition and reduction procedure. Following a 30-second test, the thermal conductivity of 
a sample was determined in less than 15 minutes: Only a comparatively small number of 
materials have been previously tested. Moreover, few materials with known thermal 
conductivities were used for validating the measurement technique. Twenty-one different 
materials, most with documented properties, were included in this investigation. As part of 
this work, the performance differences between a teflon- and glass-encapsulated thermistor 
were investigated. Kravets and Larkin used glass-encapsulated probes. Woodbury and 


Thomas used teflon-encapsulated probes. Finally, an investigation into a variation of the 
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thermistor probe (ie., the two-sensor probe) was initiated in this study. No such alternative 
to the thermistor probe has been previously reported. 

A literature review follows. Limited comparisons between the thermistor-based 
technique and other methods for measuring thermal conductivity are presented. Since 
thermistors were the only probes experimentally tested, the body of this report focuses on the 
thermistor-based measurement technique. The analytical model, the test development, and 
the results of the thermal conductivity measurements are featured. The use of the larger, 
encapsulated thermistors and the all-digital system used in this study has yet to be validated 
for measuring thermal diffusivity. For this reason, this report emphasizes the application of 


the technique for measuring thermal conductivity. 
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Literature Review 


An understanding of the other methods used to make thermal conductivity 
measurements is beneficial. Obtaining documented thermal conductivity data for calibration 
and validation purposes is critical. As is discussed later, two properties of the thermistor are 
determined by conducting calibration tests on materials with known thermal conductivities. 
Knowledge of the probe’s properties is required before the probe can be used to make 


thermal conductivity measurements of other materials. 


Other Measurement Methods 


Many methods exist to measure the thermal conductivity of liquids and solids, including 
insulation materials. Three comprehensive sources for expanded, albeit selective, 
presentations of the more common methods are references 10,11, and 12. Many technical 
papers are concerned with only one method [13,14.15,16,17]. With respect to the development 


of the thermistor probe, an understanding of four particular methods -- the guarded hot plate, 
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the flat-screen tester, the transient hot-wire, and the line source probe -- was judged to be 
important. 

A review of the methods used for measuring the thermal conductivity of insulation 
revealed that the two dominating are the guarded hot plate [18] and the flat-screen tester [19]. 
Both apparatuses are configured to promote one-dimensional conductive heat flow. The 
guarded hot plate is a steady-state measurement technique while the flat-screen tester is a 
transient measurement device. The accepted guarded hot plate standard is housed at the 
National Bureau of Standards in Gaithersburg, Maryland. This standard apparatus is reported 
to measure the apparent thermal conductivity of low-density, glass-fiber insulation with an 
uncertainty of less than 1 per cent [20]. For similar materials, the Oak Ridge National 
Laboratory flat-screen tester possesses a measurement uncertainty of + 1.7 per cent [21]. 

The probe-based technique and the two established, absolute measurement techniques 
would act in complementary roles. The guarded hot plate and the flat-screen tester would be 
used to first accurately measure the thermal conductivity of particular materials, some 
possibly being accepted as standard reference materials. These materials would be used to 
calibrate the thermistor probe. The probe would then act as a secondary standard in the field 
where cost and lack of mobility exclude the use of the two primary methods. 

The transient hot-wire and the line source designate two different types of apparatus. 
Both methods are based, however, on the same analytical model. This model consists of an 
infinitely long line source conducting to an infinite, isotropic medium. The infinite line source 
effects a constant heat flux boundary condition on the surrounding medium. In practice, the 
temperature rise of the line source or hot-wire is measured as a function of time during the 
heating process. The thermal conductivity of the surrounding test medium is determined from 
the slope of the linear plot of the temperature rise vs the natural logarithm of time [22] The 
differences between the transient hot-wire apparatus and the line source probe stem from 
their respective constructions. The transient hot-wire is constructed to promote high 
measurement accuracy while the line source probe is configured to provide in situ 


measurement capabilities. 
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For the transient hot-wire apparatus, a precision resistance wire acts as both the heating 
element and the temperature sensor [23]. This resistance wire stretches along the centerline 
of a sealed vessel. The test material, usually a liquid or gas, is poured into the vessel and 
then heated or cooled to a desired uniform temperature. Once the sample has reached 
thermal equilibrium, the hot-wire is heated, and the transient temperature measurements 
begin. Total data acquisition times range from one [24] to ten seconds [25]. The short test 
interval, along with physical modifications to the apparatus and mathematical corrections, 
helps to minimize natural convection, radiation, and stray conduction losses [26,27,28]. 
Thermal conductivity measurements of liquids using transient hot-wires are subject to overall 
uncertainties ranging from +0.5 per cent to + 1.5 per cent [10]. The transient hot-wire, which 
lacks in situ. measurement capabilities, would complement the thermistor-based 
measurement technique in the same way as the guarded hot plate and the flat-screen tester. 

For a typical line source probe, a separate resistance heating wire and a temperature 
sensor are encapsulated in a metal sheathing. The temperature sensor, either a 
thermocouple [29] or a thermistor [30], is positioned at the midpoint of the cylindrical probe. 
Because the probe is rigid and is sealed, it may be used with liquids or solids. Total data 
acquisition times may range from 30 seconds to 15 minutes. Only a relatively small window 
of this total time, however, is typically used to determine the thermal conductivity of a sample 
[31]. With respect to this window of pertinent data, Drotning and Tormey [32] report that an 
incorrect result is possible if the first of two linear regions of the characteristic plot (ie., 
temperature rise versus the natural log of time) is used. Thermal conductivities accurate to 
+4 per cent have been reported by Asher [30] for a line source probe (referred to as a 
transient needle probe in his work). Like the thermistor, a line source can be used for making 
in situ thermal conductivity measurements for insulating materials [33,34]. 

The line source probe and the thermistor probe technique each have distinct advantages 
and disadvantages. For example, the thermistor makes only relative measurements. — In 
comparison, the line source probe can make both relative and absolute measurements [30]. 


An absolute measurement technique is preferred because calibration is not required [10]. 
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Moreover, when the application requires that the line source be calibrated (to account for the 
finite heat capacity and thermal conductivity of the probe), only one reference material is 
needed [30]. The thermistor probe, on the other hand, requires at least two reference 
materials for calibration. A preferred characteristic of a small thermistor probe is its ability 
to approach a point measurement. Th? thermal conductivity of a sample as measured by a 
line source, on the other hand, is a value averaged over the length of the probe as well as 
radially. Because probe operation is based on conducting to an isothermal medium, the line 
source probe is not suited for applications in materials where a temperature gradient exists. 
Consequently, a thermistor probe is better suited to measure the thermal conductivity of 
insulation in an exterior wall or ceiling. Finally, with the thermistor-based technique, thermal 
diffusivity can also be determined from the same test. The line source probe is not known to 


have been used to measure thermal diffusivity. 


Reference Materials 


The search for appropriate reference materials (ie., materials with well-documented and 
readily repeatable thermal conductivities) began at the inception of this investigation and still 
continues. Because the thermistor probe cannot readily penetrate dense solids, the group of 
applicable reference materials is comprised of liquids, soft solids like fibrous insulation, and 
substances which can be molded around the probe. Materials with thermal conductivities less 
than that of water (0.6 W/m-°C) are sought. Suitable materials with thermal conductivities 
bracketing the range of test specimens are also essential. 

Some reference materials have been investigated by a number of researchers while 
others by only one. Using a transient hot-wire apparatus, Nieto de Castro, et al. [35] and 
Takizawa, et al. [36] measured the thermal conductivity of a variety of liquids, including 


distilled water, toluene, n-heptane, glycerin, ethyl alcohol, acetone, and methanol. A line 
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source probe was employed by Aster [30] to measure the thermal conductivity of tert butanol. 
Touloukian et al. [11] combined the work of many workers and published the most extensive 
compilation of thermal conductivity data. Many other sources document thermal conductivity 
data for a range of different materials [37,38,39,40,41,42,43,44,45]. 

Uncertainties in reported property values for reference materials exist. Differences in 
some reported values are significant. For example, of four sources listing the thermal 
conductivity of ethylene glycol [11,43,46,47], values for a temperature of 25°C range from 0.251 
[46] to 0.283 W/m-°C [43], a 13 per cent difference. Moreover, two sources [11,46] report that 
the thermal conductivity of ethylene glycol increases with temperature while the other two 
sources [43,47] report the opposite trend. With any material for which more than one source 
was identified, the author attempted to pick the best value based on the method used and the 
most recent work. For the cases based on a single source, the reported values are used but 


it is recognized that experimental error is likely. 
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Analytical Model: Thermistor Probe 


The encapsulated thermistor approximates the ideal small-volume probe. Besides 
possessing the dual capablities of heating and temperature sensing, an ideal probe would 
have a negligible mass and a spherical shape, would be electrically insulated, and the 
measured temperature would be that of its isothermal surface. The operation of the 
thermistor departs from the ideal. For example, the surface temperature of the thermistor’s 
encapsulation is not measured: instead, the measured temperature corresponds to a 
volume-averaged value for the thermistor bead. The temperature gradient within the 
thermistor must be taken into account in the analytical model. The resulting heat transfer 


problem is one with coupled heat conduction equations. 


Coupled Heat Conduction Equations 


The analytical model on which the measurement technique is based takes into account 
the time-dependent temperature distributions in both the encapsulated thermistor and the 


surrounding medium. This analytical model is mathematically expressed by equations 1 and 
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2. The heat transfer mechanism is assumed to be pure conduction. The thermistor bead and 


the encapsulation (pictured within a later section on the experimental apparatus) are 


aggregated and modelled as a perfect sphere of radius a with constant thermal conductivity, 


k, , and thermal diffusivity, «,. Tne surrounding medium is modelled as isotropic and infinite 


while having similarly constant thermophysical properties, k,, and a,,. Heat generation per unit 


volume within the thermistor, q(t) , is assumed to be spatially uniform. Throughout the 


mathematical analysis, the subscript b identifies properties of the bead thermistor while the 


subscript m denotes variables associated with the surrounding medium. 


eT oT 
ih at eS Ch) ye ey O<r<a 
re or or ky ety oly 

T re 
tile Se! 5 ua aie ae 
re? or or Sm = at 


- Four boundary conditions accompany these coupled heat conduction equations. 


symmetry: 


sak (0) = 0 t>0 
no contact resistance: 
T,(a,t) = T,,(a,t) t>0 
energy balance: 
; ree ay amet ro 
infinite medium: 
Tmloo, t) = T; t>0 
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(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Initially, the temperature throughout the thermistor and the surrounding medium is constant, 


ie., 


T,(r,0) 


T 
—t 
oO 
A 
- 
A 
® 
a 


Tr(r.0) = T,; r>a_ (8) 


In order to solve for the temperature distributions, the form of the heat generation 


function, q(t), must be determined. The approach used for making this determination follows. 


Heat Generation Function 


Three different approaches have been previously used to show that in order to achieve 
a step temperature change the heat generation per unit volume must vary linearly with the 


inverse square root of time 
219 bei a a t>0O (9) 


where T and f are referred to as the steady-state and transient constants, respectively. 
Woodbury [1] showed that in order to effect a step change on the surface of the thermistor 
bead -- an ideal condition -- the bead must dissipate power according to equation 9. Assuming 
that the function is composed of a constant plus a transient term, Valvano [5] deduced the 
form of the transient term, £ t-'/? . These two analytical approaches explained what had first 
been experimentally asceriained by Balasubramaniam and Bowman [4,6]. 

An experimental approach confirmed that the encapsulated thermistors used in this 
study generate heat in accordance with equation 9. Typical temperature response curves for 
the glass- and teflon-encapsulated thermistors dissipating heat to glycerin are given in Figs. 


1 and 2. The glass-encapsulated thermistor completes its initial transient response in slightly 
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Analytical Model: 


over a second whereas the teflon-encapsulated thermistor requires approximately 3 seconds. 
The mechanism for controlling the thermistor temperature response is described later. As 
shown in Figs. 3(a) and 4(a), once the desired temperature rise is obtained the power, here 
expressed as the voltage squared, decreases rapidly with time. (Voltage squared is directly 
proportional to the power because the resistance is constant corresponding to the controlled 
elevated temperature.) Plots of the voltage squared versus the inverse square root of time 
exhibit the expected linear characteristic for both the glass-encapsulated, Fig 3(b), and the 
teflon-encapsulated thermistors, Fig. 4(b). The y-intercept and slope from these linear plots 


are directly proportional to IT and f, respectively. 


Solution Methodology 


The procedure for solving the coupled heat conduction equations is qualitatively 
described. The final expressions for the temperature distributions in the aggregate thermistor 
bead and the surrounding medium are included. The complete solution is detailed by 
Balasubramaniam and Bowan [4] and by Valvano [5]. 

Following a simple change of variables oe = T — T,) , the two second order, partial 
differential equations (equations 1 and 2) and the four boundary conditions are transformed 
into the LaPlace domain. The resulting second-order, ordinary differential equations are 
solved, and the associated boundary conditions are applied. While in the LaPlace domain, the 
expression for the temperature distribution of both the bead and the medium arc written such 
that steady-state (I) and transient (8) constants are separated. The inversion of these 
expressions to the time domain requires contour integration. The portion of the expressions 
dependent upon the steady-state constant was first inverted by Goldenburg and Tranter [48]. 


The inversion of the expressions dependent upon the transient constant contribution is 
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Power Dissipation Characteristics of the Glass-encapsulated Thermistor when 
Approximating a Step Temperature Rise 
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Figure 4. Power Dissipation Characteristics of the Teflon-encapsulated Thermistor when 


Approximating a Step Temperature Rise 
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described by Balasubramaniam and Bowman [4]. The resulting temperature distribution 


equations for the thermistor and the surrounding medium are 
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y is the variable of integration. 


Thermal Conductivity Equation 


Because of the negative exponentials, only the first three terms in equation 10 remain 
when a steady-state condition exists (ie., t-+0co). The resulting steady-state temperature 


distribution within the thermistor bead as a function of radial position is given by 


2 k 2 
ea ELIE uke Sot ot aie 
iit) ese i 3 esl 6 (. =) (12) 


This theoretical temperature distribution must be related to the experimentally 


measured thermistor temperature, 7,. In practice, the temperature measured using the 
thermistor is an average value which corresponds to the least resistive path within the 
thermistor’s ceramic bead. To correlate the analytically determined temperature with the 
measured temperature, equation 12 is averaged over the entire volume of the sphere. 


Mathematically, 


a 2 k 
pe —-= | Ty (n) 4nr°dr = T, + Mee E + +| (13) 
(4/3)za~ 0 
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or expressed as a volume-averaged temperature rise, 


Ta?| kp 1 
aS, eae ee aun 4 
AT nace + q (14) 


where AT = AT, = AT,. The constant I (=q,,) represents the steady-state heat generation 
per unit volume. The volume is that of a sphere with radius a. The heat generation equals 
the power (given by the voltage squared divided by resistance) required to maintain the step 
temperature change within the thermistor. Making the proper substitution for the steady-state 
constant, equation 14 is rearranged to express the reciprocal medium thermal conductivity 


as the dependent variable. 


4 
— = 4na - — (15) 
k 2 k 
m Ves ) b 
Thus, a relatively simple algebraic equation for evaluating the thermal conductivity of a 
test medium results. The temperature change is an input variable while the steady-state 
voltage squared is determined by extrapolating the transient data of the measured voltage, 


squared. The two thermistor properties, a and k, , are determined via calibration. 
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Experimental Apparatus and Test Development 


The primary functions of the data acquisition and control system are to measure and 
regulate the temperature of the thermistor probe. The experimental system operates in two 
different modes. One mode of operation is conventional temperature measurement where the 
heating of the thermistor by the signal current is negligible. In the other mode, the thermistor 
is deliberately heated to a controlled temperature a few degrees above that of the test sample. 


The equipment and procedures for these two tasks are described in this section. 


Thermistor Probes 


A glass-encapsulated and a teflon-encapsulated thermistor were evaluated for 
measuring thermal conductivity. The glass-encapsulated thermistor was a Fenwal GB32P8 
standard glass probe with a nominal diameter of 2.54 mm and resistance rating of 2000 22 at 
25°C. The teflon-encapsulated thermistor was a YS! 44105 precision thermistor with a nominal 
diameter of 2.80 mm and resistance rating of 3000 Q at 25°C. Although the active bead for the 


teflon-encapsulated thermistor is larger, both thermistor beads resemble prolate spheroids 


Experimental Apparatus and Test Development 23 


as shown in Fig. 5. The encapsulations are significantly different. The glass encapsulation is 
thinner and more closely approaches the desired geometry of a sphere while the teflon 
encapsulation resembles a cylinder with a “Squared-off” end. Heat shrink tubing was used to 
seal the leads against moisture by effectively extending the encapsulations. To better seal the 
shrink tubing to the glass-encapsulated thermistor, an insulating varnish was first tried. 
Because the varnish deteriorated, however, a fine-thread was later wound at the section 
where the shrink tubing over-lapped the glass stem. A third thermistor was used for 
preliminary thermal conductivity tests and is briefly referred to later. With the exception of its 
resistance rating (1000 $2 at 25 °C), this third probe was exactly like the above described 
teflon-encapsulated thermistor. Unless otherwise specified, any reference to a 
teflon-encapsulated thermistor regards the nominal 3000 22 probe. When a reference is made 
concerning “two” thermistors, the 2000 £2 glass- and 3000 22 teflon-encapsulated thermistors 
are involved. 

Rather than rely on the manufacturer’s data, the resistance vs temperature 
characteristics of the two encapsulated thermistors were experimentally determined. A Fluke 
Platinum RTD acted as the temperature standard. The measurement uncertainty of the RTD 
is estimated to be less than +0.03° C. The Fluke instrument has a digital display resolution 
of 0.01°C and an associated stability of +0.01°C. Resistance was measured using the 
identical hardware as was employed in the thermal conductivity test. Using a constant 
temperature water bath, the probes were calibrated from 18°C to 30°C using approximately 1 
°C increments. The two encapsulated thermistors were calibrated twice -- approximately four 
months apart -- to give an indication of probe stability. 

A plot of temperature vs resistance for the glass-encapsulated probe is given in Fig. 6(a). 
The inverse relationship and the relatively high sensitivity are typical thermistor 
characteristics [49]. Three different functions, including the Steinhart-Hart equation [50], 


yielded curve fits with essentially equal closeness. The fit with the simplest form, 


T = B + D In(R) (16) 


Experimental Apparatus and Test Development 24 


payejnsdeaua-uojyjay pue (do}) payejnsdesua-ssej5 :saqoig J0}siuayy “G ainbyy 


G 


T onaw 


atu 


9 £ V c 


a Or! 


25 


Experimental Apparatus and Test Development 


(°C) 


TEMPERATURE 


1488 1552 1788 18652 228288 2158 2388 


RESISTANCE (ohms ) 


(a) 


Gigs 
nw w 
@ [a] 
Q Las] 
T TS 


~ 
n 
oO 

r 


TEMPERATURE 


16.8575 74 74 7.5 rats 72 A Sag 


NATURE LOGs OFS RESISTANCE Cln (Comms) ] 


(b) 


Figure 6. Thermistor Characteristics: (a) Temperature Versus Resistance (b) Curve Fit for 
Temperature Measurements 


Experimental Apparatus and Test Development 26 


was actually used. Figure 6(b) shows a typical result when the thermistor temperature versus 
resistance data is fitted using equation 16. The constants were evaluated using the method 
of least squares [51]. The standard error for both probe curve fits is estimated to be less than 
+ 0.04°C. From a comparison of the curve fits from the data sets collected four months apart, 
the repeatability was + 0.01 °C for the teflon-encapsulated thermistor and + 0.03 °C for the 


glass-encapsulated thermistor. 


System Hardware 


The data acquisition and control system is more readily described once the two mode 
Operation of a thermistor probe is understood. The configuration of the experimental 
apparatus and the development of the programming logic for the thermal conductivity test was 
aided by a related investigation by Fanney and Dougherty [52]. 

A thermal conductivity test begins by embedding the probe into an isothermal test 
sample. Enough time is allotted to insure that the probe-sample system achieves thermal 
equilibrium. Shortly before the actual transient measurements begin, the thermistor is used 
to measure the temperature of the sample. From this ambient temperature measurement and 
the specified temperature rise (either 2 or 3°C), the computer calculates the target thermistor 
temperature. The temperature of the probe is quickly elevated to this target value by resistive 
heating within the ceramic thermistor bead. The desired thermistor step temperature 
response is approached by using digital feedback in the form of an integral controller. The 
controlled variable is the current supplied to the probe. The thermistor voltage drop, tne 
output current, and the resulting thermistor resistance are recorded as functions of time. The 
test run is terminated after a preset time, usually 30 seconds. The computer with the 
appropriate statistical software is subsequently used to reduce the large amount of data 


collected. 
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Two different methods were used to achieve an initially isothermal test sample. The 
liquid samples and the powdered sulfur, which were contained in 125 ml glass bottles, were 
submerged in a Thermomix constant temperature water bath. For the insulating materials, a 
constant temperature chamber was used. A Dynatech R-MATIC heat flow meter thermal 
conductivity instrument doubled as the constant temperature chamber. As implied by the 
instrument’s name, the R-MATIC is conventionally used for measuring the thermal 
conductivity of insulation materials. The R-MATIC (which was used to give reference thermal 
conductivity values for the insulation materials tested using the thermistor probe) operates in 
compliance with an ASTM standard [53]. This operation, in general, involves sandwiching an 
insulation sample between isothermal “cold” and “hot” plates. Thermal conductivity is 
deduced from knowledge of the plate temperatures, the heat flux acrossed the sample, and 
the width of the sample once a steady-state heat transfer condition exists. When used as a 
constant temperature chamber, the two isothermal plates of the R-MATIC were operated at 
the same temperature. For all of the tests, the same set temperature of 25°C was used. 

The centerpiece of the experimental apparatus was a Hewlett-Packard 9836A Series 200 
computer. As is implied by Fig. 7, the 9836A computer controlled the other hardware 
components, including the thermistor. Besides !/O control, the computer 

(1) prompted the user for pertinent test parameters (eg., desired temperature 


rise, total run time, etc), 


(2) executed the integral control logic, 
(3) recorded the pertinent run variables as a function of time, and 
(4) printed and stored to a magnetic disc the collected data. 


The internal clock of the 9836A was used for the time measurements. The time measurements 
using the crystal controlled clock are accurate to one one-hundredth of a second [54]. 

The DAC Power Supply Programmer (HP 59501B) and the Precision Constant Current 
Source (HP 6186C) controlled the electrical current supplied to the thermistor. After executing 


the feedback logic, the computer sent a digital word representing the desired output current 


Experimental Apparatus and Test Development 28 


Computer 


OLLELELL LT Ld 
S © © 2 |—Power Supply Programmer 


Data Acquisition System 


Relay Box Ilencopsutatec 
Thermistor 


Figure 7. Configuration of Experimental Apparatus 


Experimental Apparatus and Test Development 23 


to the HP 59501B. The HP 59501B converted this digital word into an analog voltage signal. 
The HP 6186C, which was physically hard-wired to the DAC Power Supply Programmer, 
sensed the change in this analog voltage and updated accordingly the dc current it supplied 
to the thermistor load. 

At the start of a thermal conductivity test, the zero and full-scale current outputs from the 
remotely controlled power supply were set. A Hewlett-Packard Digital Multimeter (HP 3466A) 
was used as the amperage standard. The accuracy of the multimeter was checked against a 
lab standard and found to be within + 2 vA over the range of operation. 

When operated in the remote mode, the specified accuracy of the HP 6186C is + 0.5 per 
cent of the output current [55]. Using the HP 3466A multimeter, the remotely operated HP 
6186C power supply was tested for repeatability. In both this test and applications with the 
thermistor, the HP 6186C exhibited excellent precision, + 1A. The resolution of the remotely 
controlled current source was limited by the digital-to-analog conversion performed by the 
HP 59501B. Specifically, the interval from zero to full-scale is evenly divided into 999 steps. 
For a typical test run, a full-scale current of 3.000 mA had an associated current output 
resolution of 3 pA. 

The Hewlett-Packard Data Acquisition and Control Unit (HP 3497A) made both the initial 
thermistor resistance (ie., temperature) and transient voltage measurements. Resistance was 
measured using a 4-wire arrangement and a 10 wA signal current. This signal current caused 
negligible self-heating of the thermistor. For such resistance measurements, the resolution 
is 0.005 per cent of the reading, and the uncertainty is estimated to be less than + 0.05 per 
cent of the reading [56]. The voltage measurements were made by the resident digital 
voltmeter (DVM) of the HP 3497A. The DVM incorporates a multislope A/D converter which 
integrates the voltage signal over time. Of the three options for the measurement resolution, 
the maximum, 5 1/2 digits, was used exclusively. The software-controlled voltmeter autozero 
function was disabled for the majority of the tests before initiating the transient voltage 
measurements. When the autozero function is enabled, the DVM better compensates for any 


internally generated offset errors but requires more time to make a measurement. The 
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function was usually disabled because the offset error was believed to change minimally for 
the 30-second test and integration and 5 1/2 digits promote accuracy and resolution at the 
expense of the sampling rate. For cases where the autozero was enabled, the uncertainty of 
the transient voltage measurements is estimated to be less than +0.008 per cent of the 
reading [56]. 

The dummy resistor, housed in the relay box pictured in Fig. 7, provided a load for the 
HP 6186C power supply while it obtained its initial current setting (usually full-scale). 
Immediately before the transient voltage measurements began, the current supplied to the 
dummy resistor was switched to the thermistor. As a precautionary step, the dummy resistor 
replaced the thermistor as the external load on the power supply once the transient 


measurements terminated and the current output was zeroed. 


System Software 


Software was developed to permit testing and to aid the data reduction sequence. The 
two main programs were the thermistor resistance vs temperature calibration and the thermal 
conductivity test programs. These programs were written in HP-compatible, JEM FORTRAN77. 
IEM FORTRAN/77 code is compiled and executed within the Hewlett-Packard Pascal operating 
environment. A limitation of using the IEM FORTRAN/7 for the two tasks is that the language 
lacks direct commands for I/O controlling of the data acquisition and control equipment. To 
overcome this I/O problem, the IEM FORTRAN77 programs called subroutines contained 
within a compiled HP Pascal program module. This HP Pascal program contained special !/O 
commands that are an integral part of the Hewlett-Packard Pascal language. Finally, a short 
program, written in HP BASIC, configured the test data before it was uploaded and analyzed 
using a copyrighted software package, the HP 98820A Statistical Library [51]. An IEM 


FORTRAN77 program that runs the thermistor-based thermal conductivity test, the Pascal 
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program that permits the I/O, and the BASIC program that aids in the data reduction process 
are listed in the Appendix. The procedure for running these coupled programs accompanies 


the documented listings. 


Integral Feedback 


Integral feedback permitted controlling the temperature response of the thermistor such 
that a step change was approximated. The control was implemented using digital as opposed 
to analog technology. A block diagram of the feedback system along with the mathematical 
representation of the controller are pictured in Fig. 8. Although thermistor resistance was the 
experimental feedback variable, to aid clarity, temperature rather than resistance is used here 
to explain the feedback strategy. The difference between the desired elevated temperature, 
Ty, , and the instantaneous thermistor temperature, 7, , yielded the error, e, which was 
integrated over time. While the error was driven to zero (ie., Ty = T,), the integrated error 
approached a finite constant. The gain setting, K, scaled the integrated error to give a 
correction current, /, , of a suitable, relative magnitude. The correction current was added to 
a reference current, /,. For all the tests conducted during this investigation, the reference 
current was the full-scale current. The new current setting, /, , was supplied to the thermistor 
load, a new thermistor temperature was measured, and the process repeated. 

Two problems related to the actuator becoming saturated were overcome using 
programming logic. If the control logic yielded an updated current output that was greater 
than full-scale or that was negative, then the updated output was set to the appropriate 
boundary value. Secondly, a phenomenon known as integrator windup was avoided by 
temporarily disabling the error integration [57]. Continued integration during a saturation 
condition causes potentially large integrator build up. This build up translates into overshoots 


in the system response once the error that drove the actuator to saturation changes sign. 
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When the integration is interrupted, overshoots are minimized with the control variable 


(current, for this application) being reduced at the proper instant. 


System Speed 


System speed as used here is the frequency of measuring the thermistor’s voltage and 
updating the current output. As shown below, system speed influenced the effectiveness of 
the integral controller and the amount of data collected during the short test interval. 

An acceptable balance between maximizing system speed and measurement accuracy 
was achieved for the thermal conductivity measurement technique. During this investigation, 
system speed was improved from an initial value of slightly more than 4 Hz to a final value 
of 21.6 Hz. Measures taken to increase system speed included disabling the autozero function 
(15.6 Hz with autozero enabled) and making the thermistor voltage measurement at the lone 
external port of the HP 3497A’s DVM as opposed to using a channel on a multiplexing card. 
The loss in speed resulting from using a multiplexer card channel is attributed to additional 
internal relay switching within the HP 3497A. After these changes, the main speed limitation 
was directly linked to the integrating A/D conversion performed by the HP 3497A DVM. For 
the 5 1/2 digit voltage measurements with the autozero function disabled, the maximum 
sampling rate, according to the manufacturer’s literature, ranges from 27 to 50 Hz [56]. The 
lower bound appears to be more representative of the ceiling sampling rate for the thermistor 


application. 
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System Performance: Temperature Response 


The extent to which the experimental system performance approached the theoretical 
probe operation was investigated by using the temperature response plots from thermal 
conductivity tests. The temperature response plot was the most available indicator of system 
performance. For example, the contact resistance between the probe and the surrounding 
medium (which should theoretically be zero) and the spatial distribution of the bead heat 
generation (which should ideally be uniform) are two other performance indicators that are 
comparatively more difficult to experimentally quantify. Thus, the development and 
refinement of the experimental apparatus concentrated on reducing the initial transients as 
the probe tries to approach a step temperature change. The theoretical step response can 
only be approximated because the thermistor has a finite even if small heat capacity and, in 
addition, the power supply has a finite current limit. 

Once the basic measurement apparatus and software codes were streamlined to yield 
the fastest temperature response plots, four test parameters were investigated for their 
influence on system performance. These parameters were the system speed, the gain setting 
on the integral feedback, the power supply full-scale current setting, and the resistance rating 
of the thermistor. 

The effect of system speed was initially investigated by comparing the results from a fast 
and aslowrun. The system speed results, as shown in Fig. 9 are from a test with the nominal 
3000 Q teflon-encapsulated thermistor dissipating heat to a ethylene glycol medium. The 
thermistor completed the initial transient response in approximately seven seconds when the 
system speed was 15.6 Hz as shown in Fig. 9(a). Figure 9(b) shows that the faster system 
speed (21.6 Hz) requires approximately three seconds. Since the thermal conductivity 
equation (eqn. 15) is based on a temperature change that remains constant with time (and, 


ideally, should be a step), the data collected during the initial transients were not used when 
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plotting the characteristic data of V*vst-/?. The system speed is discussed further in 
connection with the thermal conductivity test results. 

The gain setting and the power supply full-scale current had minimal influence on the 
temperature response time while the importance of the resistance rating was second to 
system speed. The gain setting and power supply full-scale current setting influenced the 
shape of the temperature vs time plot and, if either or both were made too large, the system 
became unstable. The condition of instability was experimentally observed. Specific values 
for the gain and full-scale current settings which yielded the unstable response were not 
rigorously studied because the condition was dependent upon the combination of many test 
variables (eg., the probe, the test medium, the system speed, etc.). The effects of gain setting 
and maximum current setting are exhibited in Figs. 10 and 11, respectively. For these plots, 
the nominal 1000 22, teflon-encapsulated thermistor dissipated heat to a glycerin medium. The 
system speed was 21.6 Hz. Referring to Fig. 10, the larger gain setting produced less 
overshoot than the lower setting but also yielded more fluctuations. Referring to Fig. 11, the 
larger full-scale current gave a faster temperature rise but caused greater overshoots of the 
desired elevated temperature as compared to the smaller current setting. In both cases, there 
was no discernible effect on the settling time. Comparing the four graphs, all the transient 
temperature responses end at approximately five seconds. Since the gain and current 
settings had little effect on the time required to complete the transients, a last comparison 
was made between the two teflon-encapsulated thermistors. For a system speed of 21.6 Hz 
and a surrounding medium of glycerin, the nominal 3000 22 thermistor, as shown in Fig. 12, 
completed its transient response in approximately three seconds as compared to five seconds 
for the nominal 1000 Q thermistor. This last result led to the selection of the nominal 3000 22 


thermistor as the teflon-encapsulated probe used for the remainder of the study. 
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Calibration and Application: A General Overview 


The probe thermal conductivity calibration test is actually the same as the previously 
described measurement test. For calibration, however, equation 15 (repeated below) is 
rearranged as shown by equation 17 where C and k,, are known quantities. The two unknowns 
in equation 17, a and k,, are deduced from a plot of the resulting linear equations for each 
calibration medium. Three calibration media were used with the exception of one probe 
calibration which used two materials. Ideally, the three calibration lines (equation 17) should 
intersect at a point, thus indicating that two rather than three calibration media would have 
sufficed. In practice, however, the intersection of these three lines forms a triangle with its 
centroid being taken as the best approximation of the probe’s effective properties. Such a 
calibration plot accompanies the presentation of the measurement results for each probe in 
the next chapter. Once these effective properties are known, the application of equation 15 


is straight forward with k,, as the only unknown parameter. 


ATR 
a = Ce = ie 
mi Vss b 
ATR 
Mt A eC rere Ce a ar if (17) 
ky Km ve 


Once calibrated, each probe was used to measure the thermal conductivity of a group 
of materials. The results from the measurements were used to evaluate and, to a lesser 
degree, to validate the thermistor-based measurement technique. Understandably, probe 
evaluations were aided by comparing thermal conductivity values between those measured 
using the thermistor-based technique and a second, established method. Observed trends for 


both the glass- and teflon-encapsulated thermistors, however, provided equal insight. 
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For the probe evaluations as well as the calibration tests, materials with known k-values 
were used. The calibration materials, in addition, bracketed the thermal conductivity interval 
of eventual probe application [1]. Documented sources provided the reference thermal 
conductivity data for the liquids. For the insulation materials tested, an R-MATIC heat flow 
meter thermal conductivity instrument was used to yield reference values. 

By the time the thermal conductivity measurements documented in next chapter were 
conducted, the experimental apparatus and measurement procedure had been completely 
refined. Testing was reduced to repeating a standard procedure. For example, a typical 
thermal conductivity measurement lasted 30 seconds. The follow-up data reductions, from 
which a value of thermal conductivity resulted, required less than 15 minutes. Results and 
discussion concerning performance capabilities of the two thermistors are detailed in the next 


chapter. 
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Results and Discussion 


The overall results gave both a basis to evaluate the general measurement technique 
and the two probe types. The digital implementation of the control technique was successful. 
With some restrictions, the use of the relatively large, encapsulated thermistors was also 
validated. For the two probe types, the results were surprisingly different. The 
glass-encapsulated thermistor yielded better results when used on liquids while the 
teflon-encapsulated probe worked better on insulation materials. 

For a given test medium, multiple runs were conducted. Of the significant test variables, 
the status of the autozero function and the magnitude of the desired temperature rise were 
most often changed from run to run. A temperature rise of either 2 or 3 °C was used. The 
autozero function was usually enabled for one of the multiple runs. The standard run time 


was 30 seconds. 
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Glass-encapsulated Thermistor 


The glass-encapsulated thermistor was used to measure the thermal conductivities of 
fifteen liquids. The thermal conductivities of these liquids ranged from 0.6 (distilled water) 
down to 0.1 W/m-°C (tert butanol). Using aqueous ethylene glycol solutions, the 
glass-encapsulated probe was also tested as an interpolating instrument. The designation 
“interpolating instrument” signifies calibrating the thermistor using two reference materials 
which bracket the range of application. The calibrated probe is then used only on materials 
having compositions similar to the calibration media and thermal conductivities which fall 
within the bracket. The probe interpolates between the high and the low reference values for 
such applications. The probe was not used to measure the thermal conductivity of insulation 


materials for reasons discussed later. 


Probe Calibration 


For the glass-encapsulated thermistor, the three calibration media used were distilled 
water, glycerin, and toluene. Distilled water and toluene bracket, while glycerin divides, the 
specified thermal conductivity measurement interval. For all three materials, reliable sources 
exist which document the thermal conductivities at 25°C. The thermal conductivity test was 
repeated four times for each medium. 

Probe calibration centers on the results of two plots. From the collected data, a plot of 
V? vs t-'/ is developed and correlated with a straight line. Figure 13 shows a representative 
plot for all three of the mediums used to calibrate the glass-encapsulated thermistor. The 
y-intercept of the fitted data corresponds to the steady-state voltage squared, V2. Substituting 
for the known quantities in equation 17, a linear expression for the reciprocal bead thermal 


conductivity, 1/k,, as a function of the probe radius, a, results. Since four calibration tests 
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Results and Discussion 


were conducted on each medium, four different equations resulted. The corresponding four 
values for C (defined with equation 17) were averaged to give a single calibration curve for 
each medium. The calibration curves for the three media are plotted in Fig. 14. From the 
centroid of the resulting triangle, the probe effective radius and the reciprocal bead thermal 


conductivity were determined to be 1.52 mm and 17.6 m-°C/W, respectively. 


Probe Thermal Conductivity Measurements 


The thermal conductivity values and related test data were used to qualify probe 
performance. For the glass-encapsulated thermistor, multiple tests were conducted using 
powdered sulfur and 12 other liquids besides the three calibration media. A representative 
sample of the test data is given in Table 1. Given that an R* upper bound of 1 represents a 
perfect fit, the very high R? values from the V? vs t-‘/? plots indicate that the data sets closely 
approximated straight lines. All of the fits corresponded to R? values of 0.997 or higher. For 
the liquids tested, the value of C (equation 17) was both repeatable and independent of the 
magnitude of the temperature change. For example, the C values for identical tests with 50 
per cent aqueous ethylene glycol, table entries 13 and 14, differ by less than 0.7 per cent. For 
the same material, the C values for a 2°C and a 3°C rise are nearly identical, 19501 vs 19532. 
The results for the powdered sulfur were comparatively poor as is evident by the scatter of its 
C values. As compared to the liquids, tests on powdered sulfur yielded relatively low V2, 
values. The slower sampling rate and more accurate measurement associated with an 
enabled autozero function did not influence the V2, values in any consistent manner. Identical 
tests on glycerin (6 and 7) and toluene (10 and 11) show slightly lower values for V2, for the 
disabled condition whereas tests on acetone (2 and 3) show an opposite effect. 

The influence of system speed was further investigated by varying the sampling rate 
while testing the glass-encapsulated thermistor in glycerin. The lowest sampling rate, 2.86 


Hz, required approximately eight seconds to complete its initial transients. Therefore, in an 
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Results and Discussion 


Table 1. Glass-encapsulated Thermistor Test Data 


1.D. fie AT Autozero ve R? & 
No. Medium (°C) (°C) Function (V?) (°C/W) 
4 acetone 24.95 | Disabled 9.4696 0.9971 30478 

2 24.96 2 Disabled 6.4476 0.9988 31040 
3 24.96 2 Enabled 6.5052 0.9983 30763 
4 24.96 2 Disabled 6.4332 0.9991 31117 

5 glycerin 25,07, 2 Disabled 8.7320 0.9979 22820 
6 25.08 2 Disabled 8.7344 0.9978 22806 
7 25.08 2 Enabled 8.6952 0.9982 22906 

8 25.09 3 Disabled 12.388 0.9982 23168 

9 toluene 200 3 Disabled 8.0555 0.9992 35657 
10 25.08 2 Disabled 5.4022 0.9997 36873 
a7 25.08 2 Enabled ech al 0.9996 37186 
12 25.07 2 Disabled 5.4386 0.9995 36637 
13 50 per cent 25.02 3 Disabled 14.762 0.9992 19501 
14 ethviene glycol 25.02 3 Disabled 14.850 0.9990 19379 
15 (by mass) 25.01 2 Disabled 10.275 0.9985 19414 
16 25.02 2 Disabled 10.222 0.9978 19532 
17 n-heptane 25.08 2 Disabled 4.9435 0.9997 40292 
18 25.10 3 Disabled 7.2069 09995 39808 
19 25.10 $) Disabled 7.2144 0.9990 39770 
20 Zo. 10 2 Disabled 5.0947 0.9995 39076 
21 powdered 24.95 3 Disabled 1.5892 0.9994 181,640 
22 sulfur 24.99 2 Disabled 0.8149 0.9994 245,280 
23 25.00 2 Enabled 0.7853 0.9989 245,510 
24 25.01 2 Disabled 1.0088 0.9994 197,980 
25 methanol 252 2 Disabled 8.3417 0.9996 23845 
26 25.12 2 Enabled 8.3059 0.9995 23947 
27 252 3 Disabled 11.934 0.9992 24018 
28 2om2 3 Disabled 11.959 0.9991 23972 
29 castor oil 29.07 2 Disabled 6.7400 0.9992 29566 
30 25.08 2 Enabled 6.6856 0.9988 29791 
31 25.09 3 Disabled 9.5502 0.9994 30058 
32 25.09 2 Disabled 6.7326 0.9993 29575 
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effort to better compare characteristic plots and ultimately V2, values, the curve fits for each 
sampling rate tested were determined from the data collected after 8 seconds had elapsed. 
Because the initial temperature was nearly identical for each test (25.03°C) and a 2°C 
temperature rise was always used, determining the effect of system speed was reduced to 
comparing V2, values. As was observed in Fig. 9, a higher sampling rate yields a temperature 
response more closely approximating the theoretical step change. Despite this dependency, 
the results given in Table 2 show that the V2, values were not systematically affected by the 
sampling rate. Assuming the trend depicted in Table 2 continues to hold for sampling rates 
greater than 21.6 Hz, attaining a step temperature change is not as critical as originally 
believed. As is shown later, however, a fast system speed is important for other reasons. 
The thermal conductivities of the materials measured using the glass-encapsulated 
thermistor are listed along with those of a second documented source in Table 3. The 
calibration materials, designated by asterisks, are included to quantify the influence of having 
the calibration curves form a triangle as opposed to intersecting at a point. The magnitude 
of the largest per cent difference for these calibration and application materials is 0.7. This 
result suggests that the relative size of the calibration triangle effectively approximates a 
point. For ten of the other twelve liquids, the measurements made using the thermistor 
differed by less than 8 per cent when compared to reference values. The measured thermal 
conductivities for the aqueous ethylene glycol solutions, especially for the 50 per cent solution, 
were in good agreement with the reference values. The lack of agreement between the 
measured and reference values for methanol, was unexpected. The test data for methanol 
was comparable to data of other liquids and, in fact, the precision of its C values were the best 
of any material tested, ie. the per cent difference between the high and low C values was 0.7. 
The reason for this deviation remains unexplained. The measured thermal conductivity of 
powdered sulfur is not expected to be representative of the the true value. This belief is based 
on the significant random scatter observed in the values of C for the four test runs. A 
reference k-value for powdered sulfur was never identified. Reference 41 lists the thermal 


conductivity of solid sulfur which is 0.24 or 0.27 W/m-°C (depending on its structure, 


Results and Discussion 49 


Table 2. Steady-state Voltage Squared Values as a Function of the Sampling Rate 


Sampling Rate Ve 
(Hz) (V?) 
2.86 8.4908 
4.80 8.5264 
6.21 8.5766 
8.14 8.4986 
11.81 8.5361 
17.88 8.5046 
19.56 8.5663 
21.58 8.53109 


Probe: Glass-encapsulated Thermistor 
Medium: Glycerin 
Maximum Current: 2.700 mA 


Limits on the Reduced Data: 8s<t< 30s 
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Table 3. Glass-encapsulated Thermistor Thermal Conductivity Measurement Results 


Medium 


castor oil 
glycerin“ 
toluene” 
n-heptane 

ethyl alcohol 
methanol 
powdered sulfur 
distilled water* 
acetone 
ethylene glycol 


tert butanol 


20 per cent (by mass) 


ethylene glycol 


40 per cent 
ethylene glycol 


50 per cent 
ethylene glycol 


60 per cent 
ethylene glycol 


80 percent 
ethylene glycol 


Thermal Conductivity 


(W/m « °C) 

Thermistor Reference 

k, k, 
0.181 0.180 
0.290 0.292 
0.132 0.1311 
Oe 7 0.1228 
0.192 0.168 
0.266 0.199 
0.0158 -- 
0.611 0.6067 
O14 0.159 
0.258 0.257/0.28 
O22 0.119 
01532 0.520 
0.447 0.448 
0.418 0.418 
0.376 0.39 
0.312 O30 


“designates calibration material 


**Per Cent Diff. = 
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amorphous or polycrystalline). The observed discrepancies for methanol and powdered sulfur 
support validating a calibrated probe on as many of the applied materials as possible. 

As part of an idea which evolved from experience and was further influenced by 
Woodbury [1], the probe-based technique was evaluated for making interpolating 
measurements. Interpolating means here that the probe would only be used on materials 
having compositions and thermal conductivities similar to those of the calibration media, as 
in the case of aqueous solutions or, possibly, insulation materials. Two calibration media 
would be used and their thermal conductivities would set the extreme limits of application. 
The calibration procedure is simplified to plotting only two calibration expressions (or 
equivalently, simultaneously solving two equations with two unknowns). 

The glass-encapsulated thermistor performed well as an interpolating instrument. Using 
100 per cent ethylene glycol and distilled water as the calibration media, the probe effective 
radius and reciprocal thermal conductivity measurement gave 1.31 mm and 13.92 m-°C/W. 
Given these effective probe properties, the thermal conductivities of five different aqueous 
ethylene glycol solutions were determined. The measured results as shown in Table 4 agree 
well with the reference values. The measured thermal conductivities for 60 and 80 per cent 
solutions were in complete agreement. Figure 15 shows a curve fitted to the reference values 
and the data points from the thermistor probe measurements. Reference 43 was used 
exclusively for reference k-values of all seven materials. As discussed earlier, significant 
discrepancies exist among the reported values for the thermal conductivity of ethylene glycol. 
The accuracy of the values reported by reference 43 are used even though similar 
discrepancies are likely for the aqueous solutions. Ethylene glycol solutions were used 
because they have similar compositions and thermal conductivities that span a relatively wide 
range (0.25 to 0.6 W/m-°C). Finally, the effective properties from this second calibration differ 
from those of the first (1.31 mm vs 1.52 mm and 13.95 m-°C/W vs 18.16 m-°C/W), thus exhibiting 


the influence of the selected calibration materials on the probe parameters. 
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Table 4. Probe Application as an Interpolating Instrument 


Thermal Conductivity 


(W/m « °C) Per cent 
Medium Thermistor Reference Difference 
20 per cent (by mass) 0.534 0.520 -2.7 
ethylene glycol 
40 per cent 0.459 0.448 -3.5 
ethylene glycol 
50 per cent 0.432 0.418 -3.3 
ethylene glycol 
60 per cent 0.393 0.39 0 
ethylene glycol 
80 per cent R332 0.33 0 


ethylene glycol 
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Results and Discussion 


Teflon-encapsulated Thermistor 


The teflon-encapsulated thermistor was tested using five insulation materials, powdered 
Sulfur, and ten of the 15 liquids measured with the glass-encapsulated thermistor. The five 
insulation materials included a 100 kg/m? fibrous insulation batt, a 50 kg/m! fibrous insulation 
blanket, an alumina-silicate ceramic fiber blanket, a polystyrene sample, and a polymer foam 
slab. The results for the insulating materials and the higher viscosity liquids are emphasized. 
The majority of the results for the lower viscosity liquids are purposefully omitted because of 


poor probe performance associated with testing these materials. 


Probe Calibration 


The teflon-encapsulated probe was calibrated using glycerin, castor oil, and the ceramic 
fiber blanket. Representative plots of the V? vs t-'/? are given in Fig. 16 for the three calibration 
materials. The calibration curves of 1/k, vs a are given in Fig. 17. From the centroid of the 
resulting triangle. the effective probe radius and reciprocal thermal conductivity were 
determined to be 3.66 mm and 32.4 m-°C/W, respectively. 

The choice of the three calibration media was influenced by the findings from preliminary 
thermal conductivity tests. Using the plots of V* vs t-‘? as an indicator, the results were 
significantly better for the higher viscosity liquids and the insulation materials than for the 
lower viscosity liquids. Glycerin, castor oil, and ethylene glycol, yielded the best results 
among the liquids. Of the remaining liquids, only n-heptane and tert butanol gave fair results. 
Figure 18 illustrates the range of probe performance. The data for the polymer foam insulating 
material yields the characteristic linear curve while the data for acetone exhibits the 


randomness typical of the lower viscosity liquids. 
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Figure 18. Range of Teflon-encapsulated Probe Performance: (a) Polymer Foam Slab (b) Acetone 


Results and Discussion 


58 


The choice of the ceramic fiber blanket as the calibration raierial for the lower end of 
the thermal conductivity interval was based on a secondary factor. The reference thermal 
conductivities of all five insulation materials were measured using the R-MATIC. Therefore, 
the k-value of one material was only as well-known as that of a second. Furthermore, since 
all five materials are grouped at the lower end of 0.035 W/m-°C, the final choice was based 
on which material effected the best closure around the embedded thermistor probe. The 
ceramic fiber Blanket was chosen on this basis. 

For the teflon-encapsulated thermistor more so than the glass-encapsulated probe, the 
dependency of the measurement technique on the availability and accuracy of reference 
materials was demonstrated. The elimination of the lower viscosity liquids because of poor 
probe performance significantly reduced the number of well-documented materials. In 
particular, the elimination of distilled water and toluene, the two most extensively tested 
materials, adversely affected the final results. Left with effectively three viscous liquids and 
five insulation materials, the probe calibration was further handicapped because the 
k-values of the insulating materials were grouped at 0.035 W/m-°C. Moreover, these k-values 
were estimated to be accurate to + 5 per cent as compared to quoted uncertainties of less 
than +1 per cent for distilled water and toluene [35]. The coupled effect of the insulation 
materials having higher uncertainties and poor distribution is exhibited in Fig. 19. The 
intersection point of the calibration curves for the foam polymer and the ceramic fiber blanket 
changes from 3.19 mm and 17.4 m-°C/W to 4.01 mm and 52.3 m-°C/W as a result of a 5 per 
cent error in the measured k-value of the ceramic fiber blanket (i.e., 0.0343 to 0.0326 
W/m-°C). On a per cent basis, this 5 per cent error caused a 25 per cent increase in the probe 
effective radius and a 200 per cent increase in its reciprocal bead thermal conductivity. The 
lack of reliable reference materials obviously limits the performance of the 


calibration-dependent thermistor probe. 
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Results and Discussion 


Probe Thermal Conductivity Measurements 


Multiple tests were conducted on each material. For the insulation materials, seven 
tests were conducted on the 100 kg/m? fibrous insulation batt and the oolymer foam slab while 
four tests were conducted on the other three samples. Three tests were performed on each 
liquid and powdered sulfur. A representative sampling of the test data for the 
teflon-encapsulated probe is given in Table 5. Based on the RF? values, the linear fits again 
yield good representations for the data. The best results were for the insulating materials. 
Results for a typical lower viscosity liquid are not included in Table 5 because they are 
misleading. Specifically, related to the fact that the amount of data used to determine V2. for 
the lower viscosity liquids was comparatively small, relatively high R? values resulted. 
Referring to Table 5 for materials where the autozero function was enabled, the values of V2, 
were minimally, although consistently, higher. Although the majority of C values for the higher 
temperature rise of 3°C were higher, the effect was not systematic. As an example for the 100 
kg/m? insulation batt, the high and low C values correspond to tests with a temperature rise 
of 2°C, as shown by table entries 6 and 7. Furthermore, the scatter between the seven entries 
for this material is observed to be random. Finally, for the same material, identical test runs 
5 and 6 yielded nearly identical C values of 43838 vs 43800. This result exemplifies the high 
precision of the measurements. 

The thermal conductivities of eleven of the sixteen materials tested are presented in 
Table 6. The reference values for the insulation materials are from tests using the R-MATIC 
heat flow meter apparatus on the identical samples. Referring to the results for the three 
calibration media first, the per cent differences range from -2.8 for castor oil up to 6.1 for the 
ceramic fiber blanket. These values suggest that calibrating over the relatively large thermal 
conductivity interval and using materials with different compositions degrades probe 
measurement accuracy. Despite possibly reaching the application limits of the technique, the 


results for the majority of the tests are encouraging. The calibrated probe exhibited good 
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Medium 


100kg/m? 


insulation batt 


n-heptane 


i 


(°C) 


25.02 
25.07 
25.06 
25.04 
24.99 
24.99 
25.01 


25.06 
25.06 
25.07 


ethylene glycol 24.95 


castor oil 


polystyrene 


powdered 
sulfur 


24.96 
24.95 


25.01 
25.04 
25.05 
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24.99 
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Teflon-encapsulated Thermistor Test Data 


Autozero 
Function 


Disabled 
Disabled 
Enabled 

Disabled 
Disabled 
Disabled 
Disabled 


Disabled 
Enabled 
Disabled 


Disabled 
Disabled 
Enabled 


Disabled 
Disabled 
Disabled 


Disabled 
Disabled 
Disabled 
Disabled 


Disabled 
Disabled 
Disabled 


R2 


0.9990 
0.9992 
0.9982 
0.9993 
0.9990 
0.9990 
0.9991 


0.9859 
0.9793 
0.9825 


0.9659 
0.9785 
0.8940 


0.9967 
0.9950 
0.9858 


0.9976 
0.9979 
0.9981 
0.9982 


0.9973 
0.9964 
0.9946 
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Table 6. Teflon-encapsulated Thermistor Thermal Conductivity Results 


Thermal Conductivity 


(W/m « °C) Reference Per cent** 
Medium a Pee ie No. Difference 
castor oil* 0.185 0.180 11 -2.8 
glycerin” 0.281 0.292 30,36 3.8 
n-heptane 0.251 0.1228 aS -104 
sulfur 0.072 -- -- -- 
ethylene glycol 0.262 0.257/0.28 11/43 -1.9/7.1 
tert butanol OF153 0.119 30 -28.6 
100 kg/m? 0.0389 0.0333 R-MATIC -16.8 
insulation batt 
polymer foam slab 0.0411 0.0429 R-MATIC 42 
50 kg/m? 0.0325 0.0336 R-MATIC 3.3 
insulation blanket 
ceramic fiber blanket’ 0.0322 0.0343 R-MATIC 6.1 
polystyrene 0.0303 0.0361 R-MATIC 16.1 

ne a 

“designates calibration material 
"*Per Cent Diff. = saa x 100 
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performance for the majority of the insulating materials. The magnitudes of the per cent 
differences ranged from 3.3 to 16.8. For the only other higher viscosity liquid tested, ethylene 
glycol, the measured thermal conductivity splits the two reference values. The result for tert 
butanol was expected but the large per cent difference for n-heptane was not. In fact, the 
results for n-heptane are representative of the results for all of the other lower viscosity liquids 
tested: the measured thermal conductivity was approximately twice that of the reference 
value. The measured thermal conductivity of the powdered sulfur given in Table 6 is expected 
to be more representative of the true value than the value obtained using the 
glass-encapsulated probe (0.016 W/m-°C). 

For this study, the calibration and applied tests were essentially interchangeable since 
materials with known thermal conductivities were used in this study. Related to this 
circumstance, each material was tested before any final decision was made on which 


materials would act as the calibration media. 


Probe Comparison and Limitations 


The overall performance of each probe is directly linked to its characteristic plots of V? 
vs t-'?. Basically, the three different types of materials tested using the probes were relatively 
high viscosity fluids like castor oil, comparatively low viscosity fluids like n-heptane, and 
insulation materials like the 100 kg/m? fibrous batt.. For these three materials, a characteristic 
plot for each probe is given in Figs. 20, 21, and 22. 

For the more viscous materials like castor oil, Fig. 20, both probes showed good linearity 
after the initial transients died out. The glass-encapsulated thermistor consistently achieved 
the desired temperature change faster than the teflon-encapsulated probe. Because of its 


faster response, the glass-encapsulated probe yielded useful data sooner. 
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Characteristic Curves for Castor Oil: (a) Glass-encapsulated Thermistor (b) 
Teflon-encapsulated Thermistor 
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Figure 22. Characteristic Curves for Fibrous Insulation Batt: (a) Glass-encapsulated Thermistor 
(b) Teflon-encapsulated Thermistor 
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For the liquids with relatively low viscosities like n-heptane, the comparatively faster 
responding glass-encapsulated thermistor yielded better results. Referring to Fig. 21, the data 
for the glass-encapsulated thermistor is smoother than that of the teflon-encapsulated probe. 
The linear portion of the glass-encapsulated thermistor is also more easily distinguished as 
compared to the teflon-encapsulated probe. For example, for the glass-encapsulated 
thermistor the linear portion corresponds to those points between 0.5 and 0.9 on the x-axis 
whereas the linear portion for the teflon-encapsulated thermistor is largely dependent upon 
the researcher’s judgement. 

Using the plot of the glass-encapsulated thermistor as an example, the data graphed to 
the left of 0.50 on the x-axis depicts a time after which the probe departed from its theoretical 
operation. This departure is largely the result of the initiation of free convection. When free 
convection begins, the power dissipated by the probe increases to maintain the specified 
temperature rise. Free convection was minimized by using relatively small temperature rises. 
With materials where free convection occurs, it is important to maximize the thermistor time 
response so that the transients end and the collection of useful data begins within the minimal 
time interval. 

As shown in Fig. 22 for the 100 kg/m? fibrous insulation batt, linear plots were typical for 
data collected from tests with insulation materials. Data for the teflon-encapsulated thermistor 
are observed to maintain a linear relationship after the initial temperature transients subside. 
Furthermore, even after 30 seconds, the teflon-encapsulated curve does not exhibit the 
nonlinearity of the glass-encapsulated thermistor plot. This nonlinearity is believed to indicate 
appreciable stem losses. The signal wires embedded in the thermistor bead provide a 
significant conductance for additional heat loss. As in the case with free convection, the 
power dissipated by the probe increases to offset this heat loss. the thermal conductivities 
of the insulation materials because of a consistent peculiarity. Referring to Fig. 22, the curve 
fit for the linear portion of the curve (x greater than 0.45) extrapolated to a negative intercept 
which is, of course, not possible. The trend towards negative intercepts was also exemplified 


by the plots for the powdered sulfur. For example, values of V2, for powered sulfur ranged 
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from 0.8 to 1.6 (Table 5). Using a second insulation material, the specified temperature rise 
was varied to observe how it influenced probe performance. As shown in Fig. 23, the intercept 
is more negative for the higher temperature rise of 5°C. The reason for this peculiar tendency 
of the glass-encapsulated thermistor remains unknown. 

The glass-encapsulated thermistor was not used to measure 

A final observed trend was that the larger, teflon-encapsulated thermistor, for a given 
medium, always required more power to maintain the specified temperature rise. The 
characteristic plots of V? vs t-‘/? cannot be used for direct comparisons because the glass- 
and teflon-encapsulated thermistors have different resistance ratings (200022 vs 30002, 
respectively). The value of C, however, can be used for comparisons since it is directly 
proportional to the temperature rise divided by the steady-state power, equation 17. 
Therefore, for a given medium, a higher C value corresponds to a smaller probe power output. 
Comparing the average C values for castor oil in which both probes exhibited good 
performance, the value for the glass-encapsulated thermistor is 29748 (Table 1) versus 16293 
(Table 5) for the teflon-encapsulated probe. This representative result implies that on average 
the teflon-encapsulated thermistor required approximately twice as much power as the 
glass-encapsulated probe during a thermal conductivity measurement. 

Unfortunately, the two potential reasons investigated for the differences in probe power 
requirements suggested that the glass-encapsulated, not the teflon-encapsulated probe, 
should have the higher power rating. For example, the glass-encapsulated thermistor has the 
higher effective bead thermal conductivity of 0.057 versus 0.031 W/m-°C for the 
teflon-encapsulated probe. For a given medium, the higher conductivity probe element 
appears to conduct heat to the surface of the surrounding medium faster, and, thus, dissipate 
more energy. Secondly, the diameter of the lead wires embedded in the bead of the 
teflon-encapsulted thermistor is 0.008 in. as opposed to 0.012 in. for the glass-encasulated 
bead. Although the lead wires for the teflon-encapsulated probe are made of tinned copper 
versus a nickle-iron alloy for the glass-encapsulated thermistor, the latter thermistor is 


expected to act as the better heat sink. 
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Results and Discussion 


The noted peculiarities of each probe type remain unsolved. For example, an 
investigation into the possibility of a critical ratio for the medium thermal conductivity to the 
bead thermal conductivity (which is approximately one) for explaining the negative intercept 
of the glass-encapsulated thermistor when it is applied to insulation material has been 
unproductive. Analytical results on the effect of this ratio do not explain a negative intercept 
[4]. The simplest interpretation of the observed peculiarities is that the thermistor-based 
measurement technique has inherent limitations. Depending upon the material tested, either 
free convection or stem losses, or both, limit probe measurement performance. From 
experience, the thermistor appears to be best suited for applications in viscous liquids with 
relatively higher thermal conductivities (eg., castor oil, glycerin, etc.). The intimate thermal 
contact and higher thermal conductivities of such liquids help lessen the effect of stem losses 


while the higher viscosities either delay or even prevent the initiation of free convection. 


Summary 


The two probes complemented each other in covering the range of materials tested in 
this investigation. The glass-encapsulated probe was better for measurements in liquids 
while the teflon-encapsulated probe performed better on the insulation materials. 
Furthermore, the technique appears to be well suited for interpolating measurements after 


calibration with materials bracketing the range of interest. 
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Alternative Probe Development: The Two-Sensor 


Probe 


A thermistor is used as the active probe in the thermal conductivity measurement 
technique for many reasons. Thermistors are capable of both self heating and temperature 
sensing. The desired step temperature response is adequately approximated because the 
thermisor has a relatively low mass and high sensitivity. Encapsulated thermistors electrically 
insulate the heating element from the surrounding medium. 

Thermistors also have inherent shortcomings when used for measuring thermal 
conductivity. A thermistor has spatially nonuniform heat generation and the volume-averaged 
temperature is measured. As opposed to hollow, thermistors are solid which increases the 
thermal capacity. The analytical model also aggregates the encapsulation with the active 
bead material. These analytical assumptions approximate the operation of the physical probe. 
Because of the deviations from the analytical model, an alternative small-volume probe was 


sought. 
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Description 


The two-sensor probe is a second type of small-volume probe. The name two-sensor 
emphasizes that two components, a separate heating element and temperature sensor on a 
hollow stem, are integrated to make this probe. The heating element configuration is ideally 
a sphere; however, a short cylinder also appears acceptable. The surface temperature of the 
heating element is monitored by an attached temperature sensor (e.g., thermocouple or 
thermistor). The two-sensor probe is operated in the same mode as the self-heated 
thermistor, i.e., the surface temperature is quickly elevated and then maintained at a specified 
temperature by varying the power. With the two-sensor probe, however, the probe-medium 
interface temperature is measured directly. In comparison, the volume-averaged bead 


temperature is measured in the self-heated thermistor application. 


Theory 


The physical system of the two-sensor probe embedded in a test material is analytically 
modelled as an infinite medium with a spherical cavity. The probe corresponds to the cavity. 
The boundary condition on the infinite medium is effected by the controlled heat dissipation 
of the probe. Thus, for this probe configuration, the problem is reduced to investigating only 
the temperature profile in one material as compared to two for the self-heated thermistor 


probe application. 
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Spherical Probe 


The ideal two-sensor probe would have a spherical geometry and would be capable of 
effecting a step temperature change on its surface. For this theoretical probe, the heat 
conduction equation, the boundary condition corresponding to r > co , and the initial condition 
are the same as defined for the thermistor model, ie., equations 2, 6, and 8. The only change 


is in the boundary condition at r=a. Mathematically, 
East) +=) 1a ale t>0 (18) 
Using the Laplace transform method, the temperature distribution within the medium is 
T(t) — T, = (alr) AT, erfe[(r—a)/(2/amt) | (19) 
The heat rate as a function of time at r=a is 


CT,,(a,t 
q.(t) = —kmAg Sa km AT, 4na f a a (20) 
N 


or 
| 7a pt 


As with the self-heated thermistor probe, the characteristic plot is again the power dissipation 
as a function of the inverse square root of time. The thermal conductivity of the medium is 
dependent on the y-intercept while the slope gives its thermal diffusivity. 

As observed with the thermistor, any physical probe requires a finite time to achieve the 
target temperature change. To investigate the effect of the finite rise time, the same problem 
was solved but now for two different time-dependent boundary conditions. The boundary 
temperature was modelled epproximately as exponential rise and ramp functions. The 
solution method and results are described in reference 58. When plotted on the same graph 
with the solution for the ideal step change boundary condition, both time-dependent boundary 
condition curves are observed to converge to the linear curve of the ideal case. The 


convergence suggests that the simple mathematical expressions for the thermal conductivity 
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and the thermal diffusivity derived from the ideal case may be employed for the nonideal 


physical system. 


Cylindrical Probe 


The primary reason for exploring the potential of the cylindrical two-source probe is 
because it can be readily built. A small diameter resistance wire could be easily wound 
around a supporting rod. Accomplishing uniform wire heating over a sphere, on the other 
hand, is much more difficult. 

A cylindrically shaped two-sensor probe was also investigated. The actual probe is 
visualized as being constructed of a heating element, about one diameter long, wrapped 
around a lightweight supporting rod. As with the spherical two-sensor probe, a temperature 
sensor would be mounted on the heating element. The physical system of this cylindrical 
two-sensor probe dissipating heat to a surrounding medium has been analytically modelled 
[58]. This model is used to solve for the temperature distributions in the rod and medium 
ant an implicit finite-difference technique. 

The pending analysis should indicate if the plot of the power dissipation as a function 
of the inverse square root of time has the characteristic linear shape. If such a relationship 
is observed, the possibility exists for using the mathematical results from the ideal spherical 
geometry case along with an effective radius. The probe would be calibrated to determine its 


effective radius. 
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Application Considerations 


A rather primitive two-sensor probe was used to answer some initial application 
questions. This probe consists of a 3-mil type-E thermocouple attached to a 
teflon-encapsulated thermistor. Using a fine thread winding, as shown in Fig. 24, the position 
of the thermocouple junction is centered on the encapsulated bead. The basic test 
measurement procedure is the same as for the thermistor. For the new application, however, 
the thermistor is used only as a heating element. The attached thermocouple is now used for 
temperature sensing. 

A problem with using the thermistor, or a heating element with similar characteristics, 
was realized from the initial tests. Since thermistor resistance is temperature dependent, its 
power dissipation cannot be accurately determined from only knowing the current. 
Consequently, the thermistor voltage and the thermocouple emf were both measured. The 
overall result was a reduction in the sampling rate from 21.6 Hz to 3.9 Hz. As observed with 
the thermistor, the effectiveness of the control strategy deteriorates, and the amount of 
collected data is greatly curtailed. To limit measurements to temperature sensing, the active 
heating element of a two-sensor probe should have a very small temperature coefficient. A 
second, although less substantiated result, was that a thermocouple may lack the resolution 
and precision needed for sensing the relatively small surface temperature rise of 2°C. 
Another type of thermistor, possibly a chip type thermistor, may be a better alternative. 

High resistivity, fine diameter, and an insulating coating are the desired characteristics 
of the heating element resistance wire. Maximum current output and heating element total 
resistance set the maximum power rating of the heating element. The resolution of the 
Current output from the HP power supply is dependent upon its maximum current setting. In 
order to maximize current resolution, the largest resistance element, and thus the wire with 
the highest resistivity, is sought. Fine diameters, which are associated with high resistivities, 


also help to limit the total thermal mass of the probe. A wire with an insulating coating would 
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Alternative Probe Development: The Two-Sensor Probe 


permit closely spaced winding and use of the probe in an electrically conductive medium. 


From a review of manufacturer’s literature, suitable wire types are available. 


Project Status 


Due to time constraints, the investigation with the two-sensor probe was limited to 
analytical investigations which are not part of this thesis. The work consisted of analytically 
evaluating the cylindrical probe and conducting an analysis to determine the wire 
requirements for both heating element geometries. Once design parameters are determined, 


a prototype probe would be built. 
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Conclusions and Recommendations 


The following conclusions summarize the capabilities of the thermistor-based 


measurement technique. 


Conclusions 


An all-digital data acquisition and control system was successfully developed. The 
integral contro! strategy closely controlled the temperature response of the 
self-heated thermistor. Following a typical 30-second thermal conductivity test, the 


k-value of the test sample was available in approximately 15 minutes. 


The sampling rate did not influence the final thermal conductivity values for 
well-behaved materials. For such materials, attaining an initial step change was 
not highly important. The response time of the thermistor was important when free 


convection in liquids or stem losses in insulation materials occurred. 
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The effectiveness of the thermistor-based measurement technique is strongly 
dependent on the availability of standard reference materials. The two effective 
probe properties are determined by calibration using reference materials. Ideally, 
three calibration materials should bracket and equally divide the thermal 


conductivity interval of interest. 


The range of the thermal conductivity interval should span a wide enough interval 
to avoid similar calibration curves. As with the tightly grouped insulation 
materials, the effective probe properties from similar calibration curves are 
strongly influenced by the uncertainties in the thermal conductivities of reference 


materials. 


The glass- and teflon-encapsulated thermistor probes complemented each other 
for the group of twenty-one materials tested. The former performed better on the 
liquids with the majority of the measurements differing by less than 8 per cent 
when compared to reference values. The teflon-encapsulated thermistor 
performed better on the insulation materials. The differences between the k-value 
of the insulation materials measured using the teflon-encapsulated thermistor and 


the R-MATIC heat flow meter apparatus ranged from 3 to 17 per cent. 


Each probe exhibited unique limitations in measuring all the test materials. This 
result supports validating a given calibrated probe with materials which bracket the 


range of interest. 


The overall best results were obtained when the probe was used as an 


interpolating instrument on materials with similar compositions. 
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Recommendations 


The development of the cylindrical two-sensor probe should be continued. The results 
will indicate the -elative advantages of pursuing this probe type vs continued refinement of the 
thermistor-based technique. The advances made using the thermistor-based technique 
should enhance the development and evaluation of the two-sensor probe. 

The identification of standard reference materials, especially bracketing the range of 
insulation materials, is needed. The effectiveness of relative measurement techniques like the 
thermistor and the cylindrical two-sensor probe are dependent on the availability of standard 
reference materials. 

The use of the thermistor-based apparatus for measuring thermal diffusivity should be 
validated. A refinement of the analytical model for the thermistor’s operation may also be 
beneficial. For example, modified coupled heat conductions equations containing a heat 
generation function which accounted for the thermistor’s initial transients could possibly be 


solved. 
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Appendix A. 


Three programs are integral parts of the thermistor thermal conductivity measurement 
technique. These programs comprise the bulk of this appendix. The listings contain a large 
number of comment statements which help explain the operation of each program. Because 
of the detail given to commenting the listings, only the procedural aspects of using these three 
programs are described. 

Each program has a unique task. Related to this task, each program is written using a 
different language. The main program, THERMAC.TEXT, embodies the majority of the logic 
used to conduct a thermal conductivity test. This program is written using IEM FORTRAN/77. 
The second program, HPIO.TEXT, is actually a module which contains subroutines enabling 
1/O between the computer and the data acquisition and control equipment. This second 
program is accessed by the main program during execution of a thermal conductivity test. 
HPIO.TEXT is written using 3.0 HP Pascal. The compiled versions of each of the two programs 
-- identified by a .CODE suffix -- are the actual programs executed by the computer. The third 
program, CONVERSION, converts an ASCII data file generated during a thermal conductivity 
test into a structured BASIC BDAT equivalent file. This conversion is required in order to 
upload the raw data into a copyrighted regression analysis. CONVERSION is written using 3.0 


HP BASIC. Each of these files will be further described within the outline of the test procedure. 
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Within the Pascal environment, many softkeys are available for accessing features of the 
system. These softkeys are identified within the command line (located at the top of the CRT) 
by the capitalized letter of the key word. For example, the command line lists the key word 
“Editor” which is accessed by pressing the E on the keyboard. Some of the key words are 
included within the explanation of the test procedure which follows. These key words are set 


apart using quotations when they are first used. 


Thermal Conductivity Test Procedure 


The FORTRAN?77 ‘thermal conductivity program is executed within the HP Pascal 
environment. Once the computer has been booted, one system default should be checked and 
possibly updated. The “What” command is used to access a listing of the system defaults 
which include, among others, the defaults for the “Compiler”, the Editor, the “liBrary”, and the 
“System” volume. With respect to running a compiled thermal conductivity program like 
THERMAC.CODE, the liBrary default must be set to the IEM FORTRAN/77 library, named 
FTNLIB. This FORTRAN library must be on-line when the thermal conductivity program is first 
"“eXecuted”. Assuming an update is needed and the library is stored on the hard disc volume 
V11, the new |iBrary default should be entered as V11:FTNLIB., where the trailing period is 
required. With the liBrary default correctly set, the operator exists the list of system defaults 
by using the “Quit” option. 

Two programs are loaded into the computer's RAM memory before the thermal 
conductivity program is eXecuted. The first program is part of the IEM FORTRAN77 library. 
Using the ”P-load” softkey, the operator enters at the prompt, V11:10. (note the trailing period 
again). The second P-loaded program is the Pascal I/O module, HPIO.CODE. Assuming this 
module is also stored on the hard disc volume V11, the operator enters V11:HPIO.CODE (note 


no trailing period after the suffix). With these two programs in RAM, the thermal conductivity 
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program is initiated by eXecuteing the program, V12:THERMAC.CODE (example assumes the 
program is stored on the hard disc volume V12). 

The thermal conductivity program permits the user to set the maximum current output 
and prompts him for specific test information. With the HP 3466A multimeter connected to the 
terminals of the power supply, the full-scale and zero outputs are set »y using the appropriate 
screw adjustments on the power supply programmer. The prompts for test information are 
adequately explained within the execution of the program. Suggestions for numerical inputs, 
like the feedback gain setting, are included with a specific prompt. Following the transient test 
interval, the pertinent data - voltage squared, thermistor temperature, and time - are stored 
on an Pascal-initialized diskette. The computer will prompt the user to place the diskette in 
the left internal disc drive, #4. The data is always stored on the diskette under the filename 
VTT. Notably, the format of this data file is a special, non-ASCII type. The format is important 
because the copyrighted statistical software package used for follow-up data reduction 


requires the data to be in a structured BASIC BDAT file. 


Data Conversion 


The raw data file is converted into its ASCII equivalent and then into its structured BDAT 
format. The conversion to an ASCII format is achieved using the “Filer” of the Pascal system. 
By pressing the F-softkey, a group of commands are accessed including the “Translate” and 
“Remove” options. The user first Translates the data file by answering the two prompts with: 
V3:VTT and V3:filename.ASC, (V3: is the default volume name given to a floppy diskette dur'ng 
the Pascal formatting routine). The chosen filename of the data file can be no more than 
seven characters in length. The suffix .ASC must be included on each new filename. In 
preparation for a second thermal conductivity test. the original data file V3:VTT should be 


Removed at this time. The data disc with the ASCII files is transformed into the structured 


Appendix A. 84 


BDAT file by using the program CONVERSION. Having rebooted the computer using HP 3.0 
BASIC, the user loads and runs CONVERSION. A message within the program instructs the 
user about the placement of the source ASCII and the target BASIC-initialized floppy diskettes. 
An initial prompt permits the cataloging of either diskette. Once the user answers the prompts 
for the source and target filenames, the conversion process is completed in approximately a 
minute. The original data is now in a form which is readily uploaded for analysis using the 


copyrighted statistical software. 


Data Reduction 


A brief review of the minimal steps required to plot the characteristic curve of V? vs 


t-1/2 using the capabilities of the statistical package are outlined below. 


The program CONVERSION is first removed from memory by executing the 


command: SCRATCH A. 


kt With the diskette labelled BASIC STATISTICS placed in the right-hand internal 


drive, execute: LOAD “AUTOST: INTERNAL”, 1. 


. Answer the first few prompts accordingly. The prompt for the type of data should 


be answered with a 1 for RAW DATA. 


The prompt for the source of the data should be answered with a 2 for MASS 
STORAGE. 
- The prompt regarding the possibility of the data file being the SCRATCH FILE 


(filename: DATA) should be answered with an N. 
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The prompt for the physical location and filename of the data file is answered 
(assuming the data diskette is in the left-hand internal drive) by entering: 


filename:INTERNAL,4,1. 


An N response should be given to the prompt regarding if the data file was original 


stored by the BASIC STATISTICS program. 


The next three main prompts should be answered accordingly. The number of 
observations originally recorded is printed as part of the original test data and also 
during the execution of the program CONVERSION. Answers the follow-up prompts 


for the names of the variables with: V**2, TEMP, and TIME. 


Answer the following three Y/N prompts with Y. The structure of the data is all 
entries for the first variable followed by all entries of the second variable, and so 


on. Answer the fourth prompt in this sequence with a Y when ready. 


A peculiarity of this uploading process is that the first observation is in error. This 
observation, which corresponds to the test status at time equal to zero, may either 


be corrected or deleted using the EDIT feature of the software. 


The TRANSFORM feature of the statistical package is used to develop a new 
variable, the inverse square root of time, given the raw TIME data. Use the 
Algebraic Option (1) and the first listed equation. Store the new variable as 


Variable #4. 


Further choices of creating subfiles or deleting observations are available using 


the capabilities of the software. 


Finally, the actual statistical analysis and data plotting is initiated by accessing the 


ADVanced STATistics software. Specifically, after pressing the ADV STAT softkey, 
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two menus will appear consecutively. In response to the first list, enter 4 for 
Regression Analysis. In response to the second list, enter 3 for Polynomial 


Regression. 


Answer the prompts accordingly. The basic sequence involves specifying the 
subfile number (if any exist), and the independent and the dependent variables. 
The user requests a plot to be made (usually on the CRT), and inputs the 
parameters for the plot. Questions and results related to the regression analysis 
follow the data plotting sequence. With respect to these follow-up questions, the 
desired order of regression is 1 and the choice of a confidence coefficient is usually 


oS; 
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Cc 1 2 3 4 5 6 7 


C234567890123456 7890123456 7890123456 7890123456 7890 123456 7890123456 789012 
PROGRAM FILENAME: V12: THERMAC. TEXT 


C 
Cc 
Cc 
C THE A IN THE FILENAME THERMAC INDICATES THAT THIS VERSION OF THE 
C THERMISTOR CALIBRATION PROGRAM IS CONFIGURED TO RECORD EVERY CALL) 
C DATA OBSERVATION. THIS PROGRAM RECORDS THE INITIAL TRANSIENT 


C RESPONSE OF THE PROBE. 


C THE C IN THE FILENAME INDICATES THAT THIS VERSION OF THE 


C PROGRAM CONTAINS AN EXCESSIVE AMOUNT OF COMMENT STATEMENTS. 


C THIS PROGRAM MONITORS THE THERMISTOR VOLTAGE, VT; AND HOLOS THE 
C RESISTANCE; R_HOT, CONSTANT BY VARYING THE CURRENT. THE CURRENT 
C IS SUPPLIED BY THE COUPLED HP 6186C POWER SUPPLY AND HP 595018 


C POWER SUPPLY PROGRAMMER. 


C THIS PROGRAM IS WRITTEN USING IEM FORTRAN 773; THE INTERFACING OF THE 
C HP 9836 AND THE DATA AQUISITION SYSTEM IS ACHIEVED BY ACCESSING A 


C FILE WRITTEN IN HP PASCAL 3.0 (FILENAME: HPIO). 


C THE FOLLOWING IMPORT COMMANDS ARE USED TO ACCESS PRE-WRITTEN MODULES 
C CIE., PROGRAMS). WITHIN IEM'S FRAMEWORK, A GROUP OF MODULES COMPOSE 
C A FILE. THE MAJORITY OF THESE MODULES ARE FOUND WITHIN THE FILES ON 
C THE IEM LIBRARY DISKETTE. A LISTING OF WHICH MODULES MAKE-UP THE 

C OIFFERENT FILES ARE LISTED IN THE IEM FORTRAN MANUAL ON PAGE 4-2. 

C INFORMATION ON THE CAPABILITIES OF THE DIFFERENT CANNED MODULES 

C MAY BE INDIRECTLY DETERMINED BY READING CHAPTER 9 OF THE SAME 

C MANUAL. 

C THE VERSION OF THE IEM LIBRARY USED FOR THIS PROGRAM IS 3215. 

C A RECENTLY PURCHASED FORTRAN77 UPDATE HAS A FEW DIFFERENCES FROM 

C VERSION 3215 AND WAS NOT USED IN ANY OF THE PROGRAMS WRITTEN BY 


C BPO. 
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THE PREFIXES V11 AND V12 CORRESPOND TO THE FIRST TWO VOLUMES ON 
THE HP 91548 20 MBYTE HARD DISC (REFER TO BPD NOTEBOOK #2). 


THE MODULE V12:HPIO.CODE IS THE PASCAL PROGRAM THAT PERMITS I/O. 


v7 0 0O 02 9 


IMPORT HPIO IN V12:HPIO.CODE 
SIMPORT SYSDEVS IN V11: INTERFACE 
SIMPORT A804XDVR IN V11: INTERFACE 
$IMPORT SYSGLOBALS IN V11: INTERFACE 
SIMPORT IODECLARATIONS IN V11:I0 
$IMPORT GENERAL_1 IN V11:I0 
SIMPORT GENERAL _2 IN V11:I0 
SIMPORT HPIB_O IN V11:I0 
S$IMPORT HPIB_1 IN V11:10 
SIMPORT HPIB_2 IN V11:10 
Cc 
C THE FOLLOWING STATEMENT IS A COMPILER DIRECTIVE THAT DIRECTS 
C THE COMPILER TO NOT FLAG SECONDARY NON-ANSI FEATURES (EG.; 
C SYBOLIC NAMES WITH MORE THAN SIX CHARACTERS). 
$ANSI OFF 
PROGRAM THERMSC 
CHARACTER*81 STR_1 
CHARACTER®5 MED 
C THE FOLLOWING EXTERNAL STATEMENT ACCESSES THREE UNIQUE HP PASCAL 
C FUNCTIONS. INFORMATION ON THESE FUNCTIONS AND THE METHOD USED 
C TO ACCESS THEM MAY BE FOUND BY READING ABOUT THE EXTERNAL COMMAND 
C AND REVIEWING PAGES 9-6,9-7;9-32, AND 9-33 IN THE IEM MANUAL. THE 
C PASCAL 3.0 PROCEDURE LIBRARY MANUAL GIVES INFORMATION ON THESE 
C TIMING RELATED FUNCTIONS ON PAGES 216 TO 219. ALL THREE OF THESE 
C FUNCTIONS ARE PART OF THE MODULE SYSDEVS. 
EXTERNAL BEEP,SETSYSTIME,SYSCLOCK 
INTEGER WORD,MED1,COUNT, ILOOP,SYSCLOCK ; FLAGG 
REAL CONST1,WWORD 
COMMON /BLOCK1/ WNWORD,COUNT, TTIME,VT,R_HOT,CURRENT,AA;BB 


COMMON /BLOCK2/ ILOOP,FLAGG 
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COMMON /BLOCK3/ VT_TIME(1200,6) 


BEEP JUST MAKES THE COMPUTER BEEP. 

SETSYSTIME RESETS THE INTERNAL CLOCK OF THE 9836 TO ZERO. 

SYSCLOCK RETURNS AN INTEGER REPRESENTING THE NUMBER OF 
CENTISECONDS SINCE THE CLOCK WAS RESET. 

THE VARIABLES WORD AND WNWORD ARE THE INTEGER AND REAL VALUES OF 
SAME VARIABLE: THE DIGITAL WORD (2000 - 2999) REPRESENTING THE 
DESIRED OUTPUT FROM THE POWER SUPPLY. 

MED1 IS AN INPUT VARIABLE WHICH DENOTES THE SURROUNDING MEDIUM 
BEING TESTED. 

COUNT INDICATES THE NUMBER OF READINGS RECORDED. 

ILOOP INDICATES THE NUMBER OF UPDATES OF THE CURRENT THAT WERE MADE 
DURING THE TEST RUN. 

FLAGG IS A FLAG TO TELL WHEN THE RUN HAS BEEN TERMINATED. 

CONST1 IS THE GAIN SETTING ON THE INTEGRAL CONTROLLER. 

TTIME IS THE TIME IS SECONDS. 

VT IS THE MEASURED SELF-HEATED THERMISTOR VOLTAGE. 

R_HOT IS THE RESISTANCE OF THE SELF-HEATED THERMISTOR. 

CURRENT IS THE AMPERAGE SUPPLIED TO THE THERMISTOR IN MA. 

AA IS THE INTERCEPT CONSTANT OF THE THERMISTOR'S T VS LN(R) FIT. 

BB IS THE SLOPE CONSTANT OF THE THERMISTOR'S T VS LNC(R) LINEAR FIT. 


VTTIME IS THE ARRAY THAT STORES THE SIX PERTINENT RUN VARIABLES. 


INITIALIZING VARIABLES. 
COUNT=0 
ILOOP=0 
FLAGG=0 
ERROR_INT=0 


ERROR_INT IS THE INTEGRATED ERROR OF THE INTEGRAL CONTROLLER. 


CLEARING THE HP BUS -- 709 IS THE HP 3497A DATA ACQUISITION SYSTEM 
AND 706 IS THE HP 59501B POWER SUPPLY PROGRAMMER 


CALL MY_CLEAR(709) 
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CALL MY_CLEAR( 


SETTING THE ZERO 


706) 


POINT AND THE FULL SCALE CURRENT ON THE 


HP 59501B POWER SUPPLY PROGRAMMER 


1 WRITE(*,2) 


THE * REPRESENT THE CRT 


2 FORMAT(66(/)) 
PRINT *,' 
PRINT #,' ' 
PRINT *,' 
PRINT %,' 
PRINT *,' 
PRINT ¥,' 
PRINT *,' ' 
PRINT *,' ' 
PRINT *,'#1 == 
PRINT *,'#2 == 
PRINT ¥,'#3 == 
PRINT ®,'#4 == 
PRINT *,° ' 
WRITE(*,'C(A\)' 
READ(*,'(I1)') 
IF (I_ADJ .EQ. 
IF (I_ADJ .EQ. 
IF (I_ADJ .EQ. 


IF (I_ADJ .EQ. 


POWER SUPPLY PROGRAMMER INITIAL SET POINT CHECK’ 


USE THE HP 3466A DIGITAL MULTIMETER IN THE DC' 
CURRENT MODE' 
(NOTE: DISCONNECT ALL PROBE LEADS FROM THE DC’ 


CONSTANT CURRENT SOURCE BEFORE CHECKING SET POINTS)’ 


> ZERO ADJUST (CONTROL WORD = 2000)' 
> FULL SCALE ADJUST (CONTROL WORD = 2999)' 
> HALF SCALE CHECK (CONTROL WORD = 2500)' 


> FINISHED: .CONTINUE WITH THE CALIBRATION’ 


) 'ENTER THE NUMBER OF THE DESIRED CONDITION. ' 
I_ADJ 

1) WORD=2000 

2) WORD=2999 

3) WORD=2500 


G) GOTO 3 


MY_WRITESTRING IS A PASCAL SUBROUTINE (FOUND IN HPIO) THAT PERMITS 


A CONTROL WORD TO BE WRITTEN TO THE HP POWER SUPPLY PROGRAMMER 


CALL MY_WRITESTRING(706,WORD) 


GOTO 1 
3 WRITE(#,2) 


2000 CORRESPONDS 
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oe oe © ee © ) 


WORD=2000 
CALL MY_WRITESTRING( 706,WORD) 


CALL MY_CLEAR( 706) 


PAUSE ' meee AT THIS TIME, RECONNECT THE PROBE WIRING ‘sevens * 
WRITE(#,2) 
REMINDER . . . . THE DATA FILE THD: VTT MUST BE CONVERTED TO ITS 
ASCII EQUIVALENT BEFORE RUNNING ANOTHER CALIBRATION CASE 
PRINT *, ° BEFORE CONTINUING ; DID YOU RENAME THE LAST DATA FILE TO ITS' 
PRINT *, ' ASCII EQUIVALENT?’ 
PRIN Tesi sa 
PRINT *, ° IF NOT, AT THIS TIME TERMINATE THIS RUN AND GO TO THE FILER AND:' 
PRIN Tet 5 tenes 
PRINT *, ' (1) TRANSLATE THE FILE THD: VTT TO A READILY IDENTIFIABLE 


ASCII EQUIVALENT’ 


PRINT *, ° CIE. WATER_1¢.ASC -- REMEMBER THE .ASC).' 
PRINT sts.00 9 
PRINT *, ° (2) REMOVE THE FILE THD: VTT FROM THE DATA DISC; THIS 


IS ALWAYS THE NAME' 

PRINT *, ' OF THE IEM FORTRAN77 GENERATED DATA FILE.’ 
PRINT *®, ‘* ° 

PRIN TOC sec 

PRINT *, ‘* ' 

PAUSE 


WRITE(#,2) 


BEGIN THE CALIBRATION PROCEDURE 


COLLECT USER INPUT DATA 


PRINT *,'#1 ==> DISTILLED WATER’ 
PRINT *,'#2 ==> AIR‘ 
PRINT *,'#3 ==> CASTOR OIL' 


PRINT *,'#4 ==> GLYCERIN’ 
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PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 


PRINT *,' 


#,'°R5 
%,'R6 
%,'S7 
*,'H8 


Hy RO 


%,'°#10 
811 
#,'°R12 
M,'R13 


,' 1G 


==> 


==> 


==> ETHYL ALCHOL' 

==> TOLUENE’ 

==> 100 ETHYLENE GLYCOL" 
==> TERT BUTANOL‘ 


==> ACETONE'* 


SULFUR‘ 
N-HEPTANE' 
METHANOL ' 
ARMAFLEX 


NBS INSULATION BATT' 


WRITE(#,'CA\)')' ENTER THE NUMBER CORRESPONDING TO THE 


+MEDIUM. 


READ(#, '(I2)* )MED1 


PRINT *, 


PRINT *®, 


WRITE(6,10) 


C UNIT 6 IS THE HP PRINTER. 


10 FORMATC'TEST MEDIUM: ‘',\) 


LE 


IF 


age 


LE 


IF 


IF 


IF 


IF 


Mg 


r= 


IF 


d= 


1g 


IF 


15 FORMAT (A) 


(MED 1 


(MED 1 


(MED1 


(MED 1 


(MED 1 


(MED1 


(MED 1 


(MED1 


(MED1 


(MED1 


(MED1 


(MED1 


(MED1 


(MED 1 
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-EQ. 
-EQ. 
-EQ. 
-EQ. 
-EQ. 
-EQ. 
-EQ. 
-EQ. 
-EQ. 
-EQ. 
-EQ. 
-EQ. 
-EQ. 


-EQ. 


1) 


2) 


3) 


4) 


5) 


6) 


7) 


8) 


9) 


WRITE(6,15)'DISTILLED WATER’ 
WRITE(6;15)*AIR' 
WRITE(6,15)'CASTOR OIL' 
WRITE(6,15) ‘GLYCERIN’ 
WRITE(6,15)'ETHYL ALCHOL' 
WRITE(6,15)* TOLUENE’ 
WRITE(6,15)'100 ETHYLENE GLYCOL' 
WRITE(6,15)' TERT BUTANOL’ 


WRITE(6,15) ‘ACETONE ' 


10)WRITE(6,15) ‘SULFUR’ 


11)WRITE(6,15) 'N-HEPTANE' 


12)WRITE(6;15) ‘METHANOL ' 


13)WRITEC6,15) ‘ARMAFLEX' 


14¢)WRITE(6,15)'NBS INSULATION BATT' 
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WRITE(#,2) 
PRINT *, ' ' 
PRINT *, ‘ ' 
PRINT *, ‘OPTIONS: #1, #3, #7, #9, #14' 
WRITE(*,'(A)') ENTER THE THERMISTOR ID NUMBER. ' 
READ(™,'(I2)")NTH 
WRITE(6,16)NTH 
16 FORMAT('THERMISTOR ID NUMBER: ‘,I2) 
C 
C THE FOLLOWING ARE THE CONSTANTS FOR A LINEAR CURVE FIT OF THE 
C THERMISTOR'S TEMPERATURE AS A FUNCTION OF ITS NATURAL LOG OF 
C RESISTANCE. 
IF (NTH .EQ. 1) THEN 
AA=206.708 
BB=-22.6781 
END IF 
IF (NTH .EQ. 3) THEN 
AA=206.47445 
BB=-22.6413 
END IF 
IF (NTH .EQ. 7) THEN 
AA=211.9974 
BB=-27.09895 
END IF 
IF (NTH .EQ. 9) THEN 
AA=212.3583 
BB=-25.14037 
ENO IF 
IF (NTH .EQ. 14) THEN 
AA=211.58059 
BB=-26.93877 
END IF 


PRINT 1% suecen. 
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PRINT #, ° * 


THE FEEDBACK GAIN IS SOMETHING THAT YOU WILL HAVE TO PLAY WITH. THE 
GAIN'S MAGNITUDE DEPENDS ON THE RELATIVE MAGNITUDE OF THE INTEGRATED 
ERROR -- IN THIS CASE, THE ERROR IS THE DIFFERENCE BETWEEN THE ACTUAL 


THERMISTOR RESISTANCE AND THE DESIRED THERMISTOR RESISTANCE. 


wetset NOTE: THE SYSTEM HAS BEEN OBSERVED TO GO UNSTABLE WHEN TOO 


HIGH OF A GAIN IS CHOSEN. 


WRITE(*,'C(A)') ‘ENTER THE DESIRED FEEDBACK CONTROLLER GAIN 
+(5.00 IS A GOOD FIRST GUESS). ' 

READ(*, '(F5.2)')CONST1 

WRITE(6,17)CONST1 

FORMAT('FEEDBACK GAIN SETTING: ',F5.2) 

PRINT *,' ' 

PRINT ™,' ' 

WRITE(#,'C(A)') ENTER RUN TIME (MINUTES) -- 0.5 IS A GOOD CHOICE. ' 
READ(*,%)RUN_TIME 

PRINT *, ' ' 

PRINT *, ‘ ' 

RUN_TIME=RUN_TIME*60.0 

WRITE(6,20)RUN_TIME 


FORMAT('RUN TIME: ',F5.2;1X;,'SEC') 


C THE FOLLOWING LINES ARE INPUTS REGARDING THE DESIRED OPERATION 


C PARAMETERS OF HP 6186C DC CURRENT SOURCE. 


WRITE(™,'°CA)') ‘ENTER THE FULL SCALE CURRENT, MA (EX. 5.000). ° 
READ (#, #) RANGE 

PRINT *, ' 

PRIN sma 

WRITE(*,'(A)') "ENTER THE INITIAL CURRENT, MA (EX. 5.000). ° 
WRITE(™,'(A)')' CUSUALLY WANT THE FULL SCALE CURRENT AT T=0).' 


READ (*, *) CONTROL 
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PRINT *, ° ° 
PRINT *, ° ° 
WRITE(6, 30) CONTROL 
30 FORMATC'INITIAL CURRENT: ‘',F6.3,1X»'MA') 
WRITE(6,31)RANGE 
31 FORMATC "FULL SCALE CURRENT: ‘»F6.351Xs‘MA') 
CCONTROL= (CONTROL /RANGE )¥999.0 
ICONTROL=INT(CCONTROL+0.5) 
THE CONTROL WORD RECOGNIZED BY THE HP 59501B MUST BE 
BETWEEN 2000 AND 2999. THUS, A GIVEN RANGE CHOICE (THE SPAN) 


IS DIVIDED INTO 999 INCREMENTS. 


MAKING SURE THE FINITE LIMITS OF THE DIGITAL WORD ARE NOT 
EXCEEDED. 

IF CICONTROL .GT. 999) ICONTROL=999 

IF CICONTROL .LT. 0) ICONTROL=0 

WORD=2000+ICONTROL 

WAORD=WORD*1 . 0 

CURRENT=CONTROL/1000.0 


CURRENT IS IN TERMS OF AMPS. 


THE VARIABLE ICONTR IS USED LATER IN THE CONTROL LOOP. 


ICONTR=ICONTROL 


WRITE(*,'(A)') ENTER THE DESIRED DELTA OHMS (WITH RESPECT TO RECORDING 
DATA).' 

READ(#,*)DELTA_OHMS 

PRINT *, ' ' 

PRINT *, ' ' 

WRITE(™,'(A)')'ENTER THE DESIRED DELTA T -- SUGGEST 2.00 C.' 
READ(™,*)DELTA_T 

PRINT *, ' ' 

PRINT ¥, ' ' 


WRITE(6340)DELTA_OHMS 
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GO FORMATC'DELTA OHMS: ',F4.1,1X, ‘OHMS') 
WRITE(6,50)DELTA_T 


50 FORMATC'DELTA T: ‘',F2.0,1Xs5'C') 


IEM FORTRAN77 VARIABLES ARE NOT EXACTLY IN THE SAME FORMAT AS THOSE 


USED IN PASCAL. SPECIFICALLY, THE FIRST CHARACTER OF A PASCAL STRING 


MUST BE THE NUMBER OF CHARACTERS WITHIN THE STRING. FOR EXAMPLE, THE 


STRING VAl HAS THREE CHARACTERS. BUT, SINCE THE V IS IN POSITION 
2» A LEADING BLANK HAS BEEN ALLOCATED BY THE NEXT STATEMENT. THE 
BLANK IS FILLED BY THE NUMBER 3 BY THE STR_1(1:1) ASSIGNMENT 

IE.» 4 - 1 = 3). ALL OF THE STRINGS THAT ARE EXPORTED TO A PASCAL 
MODULE -- LIKE THE COMMANDS TO THE HP 3497A AND THE HP 59501B -~- 
MUST BE BUILT USING THESE SAME STEPS. THE STRING STR_1 IS 


ALWAYS USED AS THE OUTPUT STRING. 


SET THE AUTOZERO FUNCTION OF THE HP 3497A TO ON. 
STR_1(2: )="VA1' 
STR_1(1:1)=CHAR(LEN(STR_1)-1) 
CALL MY_WRITESTRINGLN(709,STR_1) 
MAKE THE THERMISTOR INITIAL RESISTANCE (=> TEMPERATURE) MEASUREMENT 
FIRST TURN ON THE HP 3497A'S KNOWN CURRENT SOURCE. 
NOTE: THE VOLTAGE CARD IN SLOT 1 (STARTING FROM 0) HAS BEEN 
CONFIGURED TO MAKE G-WIRE RESISTANCE MEASUREMENTS. 
STR_1(2: )='VC1' 
STR_1(1:1)=CHAR(LEN(STR_1)-1) 
CALL MY_WRITESTRINGLN(709,STR_1) 
PAUSE ' WAIT AT LEAST 2 SECONDS BEFORE CONTINUING 
+ ‘ 
PRINT *, ' ' 


PRINT) 0’) aoe & 


C CLOSING THE CHANNELS OF INTEREST 


STR_1(02: )='AC20,30' 
STR_1(1:1)=CHAR(LEN(STR_1)-1) 


CALL MY_WRITESTRINGLN(709;STR_1) 
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CALL MY_READNUMBER(709,V_COLD) 
THE ABOVE MAKES THE DESIRED VOLTAGE MEASUREMENT. REFER TO THE 
LISTING OF THE PASCAL PROGRAM HPIO.TEXT WHICH PERMITS THIS 1/0 
AND CONTAINS THE SUBROUTINES MY_WRITESTRINGLN AND MY_READNUMBER. 
R_COLD=V_COLD/1E-5 
R_COLD IS THE THERMISTOR'S NONSELF-HEATED RESISTANCE IS OHMS. 
WRITE (#,60)R_COLD 
60 FORMAT('R_COLD =',F7.251X, OHMS’) 
WRITE(6,70)R_COLD 
70 FORMAT('COLD RESIST: ‘',F7.2,1X)'OHMS') 
TURNING OFF THE 10 MICROAMP CURRENT SOURCE OF THE HP 34974. 
STR_1(2: )='VCO" 
STR_1(1:1)=CHAR(LEN(STR_1)-1) 
CALL MY_WRITESTRINGLN(709,STR_1) 
PERFORMING AN ANALOG RESET FOR THE HP 3497A. 
STR_1(2: )='AR' 
STR_1(1:1)=CHAR(LEN(STR_1)-1) 
CALL MY_WRITESTRINGLN(709,STR_1) 
THE FOLLOWING COMMAND PUTS THE 1000 OHM DUMMY RESISTOR IN~LINE. 
STR_1(2: )="DCG)1" 
STR_1(1:1)=CHAR(LEN(STR_1)-1) 
CALL MY_WRITESTRINGLN(709,STR_1) 
CALCULATING THE INITIAL TEMPERATURE OF THE THERMISTOR 
T_COLD=AA+(BB*¥ALOG(R_COLD)) 
DETERMINING THE DESIRED ELEVATED TEMPERATURE FOR THE THERMISTOR. 
T_HOT=T_COLD+DELTA_T 
BACKING OUT THE THERMISTOR RESISTANCE CORRESPONDING TO THIS 
DESIRED ELEVATED TEMPERATURE. 
RT_FINAL=EXP(-((AA-T_HOT)/BB) ) 
WRITE(*,80)T_COLD 
80 FORMAT('T_COLD = ',F5.2,1X,'C') 
WRITE(#,90)T_HOT 
90 FORMAT('T_HOT = ',F5.2,1X,'C") 


WRITE(6,100)T_COLD 
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100 FORMAT('COLD TEMP: ‘,F5.2,1X,'C') 
WRITE(6,110)T_HOT 
110 FORMAT('HOT TEMP: ‘,F5.2,1X,'C') 
WRITE(*,120)RT_FINAL 
120 FORMAT('RT_FINAL = ',F7.251X, OHMS’) 
WRITE(6,130)RT_FINAL 
130 FORMAT('DESIRED RESIST: ',F7.2,1X,'OHMS') 
C BEGIN THE DATA COLLECTION ROUTINE 
PRINT *, °° 
PRINT *, ' ' 
C THE FOLLOWING LINE PUTS THE DESIRED INITIAL CURRENT ACROSSED 
C THE 1000 OHM DUMMY RESISTOR -- THIS ALLOWS THE CURRENT SOURCE 
C TO GET READY. 
CALL MY_WRITESTRING( 706 ,WORD) 
CALL BEEP 
PAUSE ‘READY ‘' 
C REFER TO THE HP 3497A QUICK REFERENCE SHEET FOR INFORMATION ON 
C THE FOLLOWING COMMANDS. 
STR_1(2: )='VAO' 
STR_1(1:1)=CHAR(LEN(STR_1)-1) 
CALL MY_WRITESTRINGLN(709;STR_1) 
STR_1(2: )='VRS5' 
STR_1(1:1)=CHAR(LEN(STR_1)-1) 
CALL MY_WRITESTRINGLN(709,STR_1) 
STR_1(2: )='VDS5' 
STR_1(1:1)=CHAR(LEN(STR_1)-1) 
CALL MY_WRITESTRINGLN(709,STR_1) 
STR_1(2: )=*SDO' 
STR_1(1:1)=CHAR(LEN(STR_1)-1) 
CALL MY_WRITESTRINGLN(709,STR_1) 
C THE FOLLOWING COMMAND BRINGS THE THERMISTOR IN-LINE AND TAKES 
C THE DUMMY 1000 OHM RESISTOR OUT-OF-LINE. THE INITIAL CURRENT 
C IS NOW SEEN BY THE THERMISTOR. 


STR_1(2: )='DCG,0" 
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STR_1(1:1)=CHAR(LEN(STR_1)-1) 
CALL MY_WRITESTRINGLN(709,STR_1) 
C SETTING THE SYSTEM RUN CLOCK TO ZERO 
CALL SETSYSTIME(0) 
C INNER LOOP: MEASURE VT AND CALCULATE R_HOT 
C MY_READNUMBER IS A PASCAL SUBROUTINE THAT MEASURES AND RETURNS 
C THE THERMISTOR VOLTAGE. 
135 CALL MY_READNUMBER(709,VT) 
C CHECKING THE CLOCK AND COUNTING THE NUMBER OF VOLTAGE MEASUREMENTS 
C THAT ARE BEING MADE DURING THE RUN. 
TTIME=SYSCLOCK()/100. 
ILOOP=ILOOP+1 


R_HOT=VT /CURRENT 


Cc INTEGRAL CONTROL STRATEGY 


C THE DIFFERENCE IN THE THERMISTOR'S RESISTANCE ACTS AS THE FEEDBACK 
C ERROR. RT_FINAL IS THE DESIRED RESISTANCE SETTING. 
ERROR_R=(R_HOT-RT_FINAL) 
C THE NEXT LINE IS MOST APPLICABLE DURING THE INITIAL START-UP. IT 
C HELPS THE SYSTEM TO NOT HAVE AS GREAT AN OVERSHOOT OF THE DESIRED 
C RESISTANCE. YOU WANT THE CURRENT TO GO DOWN AS SOON AS THE DESIRED 
C RESISTANCE IS REACHED$; YOU DON'T WANT A LAG CAUSED BY AN INTEGRATION 
C OF A POSITIVE FEEDBACK ERROR. THE CONTROLS PEOPLE REFER TO THIS 
C PROBLEM AS INTEGRATOR WINDUP. 
IF (ICONTR .GE. 999 .AND. ERROR_R .GT. 0.0) GOTO 200 
IF (ICONTR .LE. 0 .AND. ERROR_R .LT. 0.0) GOTO 200 
C INTEGRATING THE FEEDBACK ERROR 
ERROR_INT=ERROR_INT+ERROR_R 
C ICONTROL IS A REFERENCE STATE CORRESPONDING TO THE INITIAL 
C CONTROL WORD, USUALLY 999 -- THE LEADING 2 IS ADDED LATER. 
200 ICONTR=INT( ICONTROL+(CONST1¥ERROR_INT+0.5)) 
C THE FOLLOWING LINES MAKES SURE THAT THE MAX VALUE OF 999 AND 


C THE MIN VALUE OF 0 ARE NOT EXCEEDED. 
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IF CICONTR .GT. 999) ICONTR=999 
IF CICONTR .LT. 0) ICONTR=0 
WORD=2000+ICONTR 
WWORD=WORD*1.0 
C DETERMINING THE CURRENT BEING SENT TO THE THERMISTOR. RECALL; 
C THE CURRENT IS KNOWN; IT IS NOT MEASURED; YOU ARE CONTROLLING IT. 
CURRENT=ICONTR¥(RANGE/999.)/1000. 
CALL MY_WRITESTRING(706,WORD) 
C DETERMINING IF THE DATA FOR THIS OBSERVATION SHOULD BE STORED. 
DELTA=ABS(RT_FINAL~R_HOT) 
COUNT=COUNT+1 
IF (COUNT .GT. 1200) CALL TERMINATE 
VT_TIMECCOUNT,1)=1/(( 1000*TTIME) ¥¥.5) 
VT_TIME (COUNT » 2) =VT##2 
VT_TIMECCOUNT , 3)=R_HOT 
VT_TIMECCOUNT »4)=1000. ®CURRENT 
VT_TIME (COUNT »5)=WWORD 
VT_TIMECCOUNT ,6)=TTIME 
IF (TTIME .GT. RUN_TIME) CALL TERMINATE 
C IF THE SUBROUTINE TERMINATE HAS BEEN CALLED, THEN FLAGG WILL BE 
C SET TO 1, AND THE PROGRAM WILL STOP TAKING DATA. 
IF (FLAGG .EQ. 0) GOTO 135 
CRIN eee 
PRINT *®, ' ' 
PRINT ®, ‘* ' 
CALL BEEP 
PRINT *, ' St 3 36 3 3 FINISHED 438 


END 


SUBROUTINE TERMINATE 
C THIS SUBROUTINE IS USED TO STOP THE TEST RUN AND GIVE THE OPERATOR 
C THE OPTION OF PRINTING AND/OR STORING THE DATA. 

INTEGER COUNT, ILOOP,1,ANSWER; FLAGG,WORD 


CHARACTER*®81 STR_1 
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COMMON /BLOCK1/ WNORD, COUNT, TTIME,VT,R_HOTsCURRENT,AA,BB 
COMMON /BLOCK2/ ILOOP,FLAGG 
COMMON /BLOCK3/ VT_TIME(1200,6) 
TURNING THE CURRENT SUPPLY TO ZERO OUTPUT. 
WORD=2000 
CALL MY_WRITESTRING(706,WORD) 
THE FOLLOWING DISCONNECTS THE HP 3497A DVM FROM THE CIRCUIT 
THIS IS NEEDED SINCE A SINCE A LONE VOLTAGE MEASUREMENT IS 
BEING MADE AT THE PORT OF THE DVM. 
STR_1(2: )='DOG,1' 
STR_1(1:1)=CHAR(LEN(STR_1)-1) 
CALL MY_WRITESTRINGLN(709,STR_1) 
TAKING THE THERMISTOR OUT-OF-LINE AND PUTTING THE 1000 OHM 
RESISTOR BACK IN-LINE. 
STR_1(2: )='D04,0' 
STR_1(1:1)=CHAR(LEN(STR_1)-1) 
CALL MY_WRITESTRINGLN(709,STR_1) 
PRINT *,' ' 
WRITE (*,220) TTIME 
220 FORMAT('TIME = ',F5.2,' SEC’) 
CPS=ILOOP/TTIME 
PRINT *,° ' 
WRITE(6,230) CPS 
230 FORMAT (‘READINGS PER SECOND: ‘',F5.2) 
PRINT *,' ' 
PRINT *,° ' 
CALL MY_CLEAR(709) 
CALL MY_CLEAR(706) 
WRITE(6,300)TTIME 
300 FORMATC'ELAPSED TIME: ',F7.2,1X,'SEC') 
WRITE(6,310)COUNT 
310 FORMAT('VALUES RECORDED: ',I4) 


WRITE (6,%)* ' 
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C OPPORTUNITY TO PRINT THE PERTINENT DATA 


c 


c 


Cc 


C 


320 


330 


350 
400 


500 


510 


WRITEC#,‘CA\)')'DO YOU WISH TO PRINT THE DATA (1 ==> Y OR 


+2 ==> N)? : 


READ (#,#) ANSWER 
IF CANSWER .EQ. 2) GOTO 500 
WRITE(6,320) 


FORMAT( '1/CTIME##.5)' 6X, 'CVT_HOT)*#2°,6X,‘RT_HOT',6X, ‘CURRENT’, 


+6X5*CONTROL' 6X5 'TIME') 


WRITE(6,330) 


FORMAT('1/(CMSEC##.5)' 9X, 'V¥H2'411X5 OHMS’ 5 9X, "MA's 10X, 


+'WORD' ,8X, ‘SEC’ ) 


WRITE (6,%)° ' 
DO G00, I=1, COUNT 


WRITE(6,350)VT_TIMEC(I,1),VT_TIME(I,2),VT_TIME(I,3), 


+VT_TIMECI54),VT_TIME(I,5),VT_TIME(1I,6) 


FORMAT (2X5F6.5;9XyF8.4,6X5F8.2,7X,F5.3538X5F5.0,;6XsF6.3) 
CONTINUE 

PRAN Te sone 

PRIN TOG; ioe 

WRITE(#,'CA)*)'DO YOU WISH STORE THE DATA ON A DISC’ 
WRITEC#,*CA\)')'C1 ==> Y, ==> ND? ' 

READ( *,*) ANSWER 


IF CANSWER .EQ. 2) GOTO 600 


STORING THE DATA TO THE DATA DISK IN DRIVE #4. 


WRITE (*,510) 

FORMAT (60(/)) 

PRINT *, ‘BEFORE PROCEEDING YOU MUST BE SURE THAT THE DATA FILE’ 
PRINT *, ‘NAMED THD: VTT HAS BEEN FIRST CONVERTED TO ITS ASCII' 
PRINT *, ‘EQUIVALENT VIA THE TRANSLATE OPTION WITHIN THE FILER' 
PRINT *, ‘AND THEN REMOVED USING THE SAME FILER. ' 


PRINT ¥, ' ' 
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PRINT *, ° ' 
PAUSE 
WRITE (#,510) 
PAUSE ' mew ~=PLACE THE DATA DISKETTE IN DRIVE #4 _ sstsesee ° 
WRITE(#,511) 
511 FORMAT(66(/)) 
PRINT *, “WRITING ... .' 
OPEN(CUNIT=5,FILE="THD: VTT',STATUS='NEW' ) 
C THE DATA DISC MUST HAVE A VOLUME LABEL OF THD IN ORDER FOR THE 
C STORING PROCESS TO BE SUCCESSFUL. THD STANDS FOR THERMISTOR DATA. 
NEND3=3*COUNT 
C NEND3Z IS THE TOTAL NUMBER OF VALUES TO BE STORED. 
WRITE(5,513)NEND3 
513 FORMATCI5) 
C STORING THE THE SQUARED THERMISTOR VOLTAGES. 
DO 520, I=1,COUNT 
WRITE(5,515)VT_TIME(I,2) 
515 FORMAT(F8.4%) 
520 CONTINUE 
C STORING THE SELF-HEATED TEMPERATURES OF THE THERMISTOR. 
DO 530, I=1,COUNT 
TEMP=AA+BB¥ALOG(VT_TIME(I,;3)) 
WRITE(5,525) TEMP 
525 FORMAT(F6.2) 
530 CONTINUE 
C STORING THE TIMES (IN SECONDS) WHEN THE MEASUREMENT WERE MADE. 
DO 540; I=1,COUNT 
WRITE(5,535)VT_TIME(1,6) 
535 FORMAT(F6.2) 
540 CONTINUE 
600 FLAGG=1 
RETURN 


END 
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MODULE HPIO}3 


{ 


THIS PROGRAM IS WRITTEN USING HP PASCAL 3.0. THE PROGRAM CONTAINS 
THE SUBROUTINES THAT ALLOW THE IEM FORTRAN77 PROGRAMS TO INTERFACE 


WITH THE HP DATA ACQUISITION AND CONTROL EQUIPMENT. HP PASCAL HAS 


SPECIAL COMMANDS THAT PERMIT THE I/O. } 


SSEARCH 'V11:10.',*V11: INTERFACE. '$ 


{ 


THE DOLLAR SIGNS INDICATE THAT THE SEARCH COMMAND IS A 
COMPLIER DIRECTIVE.} 


IMPORT GENERAL_2;HPIB_2,IODECLARATIONS; 


{ 


THIS PASCAL PROGRAM UTILIZES MODULES FOUND ON THE FILES IO AND 
INTERFACE. THESE FILES, IN TURN, ARE PART OF THE IEM LIBRARY 
(VERSION 3125). REFER TO PAGE 4-2 OF THE IEM FORTRAN77 MANUAL 
FOR MORE INFORMATION ON WHICH MODULES MAKE UP THE TWO FILES. 
UPON REVIEWING THE AVAILABLE HP PASCAL LIBRARY MODULES, I FIND 
THAT THE SAME MODULES EXIST. I NEVER TRIED TO VERIFY THAT THE 
IEM AND HP MODULES ARE THE SAME. REFER TO PAGES 6;7, AND THE 
CHAPTER ENTITLED INTRODUCTION TO I/O IN THE PASCAL 3.0 


PROCEDURE LIBRARY MANUAL FOR MORE INFORMATION. } 


EXPORT 


( 


THE EXPORT COMMAND STATES THAT THE FOLLOWING TYPES AND PROCEDURES 


MAY BE IMPORTED BY OTHER PROGRAMS (LIKE THERMSC).} 
TYPE STRING80 = STRING[80]; 
PROCEDURE MY_WRITESTRINGLN(DEV: INTEGERS STR: STRING8O) 5 
PROCEDURE MY_TRIGGERV(DEV: INTEGER) $ 
PROCEDURE MY_WRITESTRING(DEV: INTEGER; VAR A: INTEGER) 5 
PROCEDURE MY_READNUMBER(DEV: INTEGER;VAR A:REAL)3 


PROCEDURE MY_CLEAR(DEV: INTEGER) § 


THE MY_ PREFIX IS OFTEN USED TO MAKE A COMMAND INTO A SUBROUTINE 


CIE., A PROCEDURE). THE PROCEDURES CAN BE NAMED ANYTHING 
THAT HELPS THE PROGRAMMER RECOGNIZE WHAT IT DOES. THE DEV 
VARIABLE STANDS FOR THE HP ADDRESS CODE. FOR EXAMPLE, THE 
HP 3497A IS 709 WHILE THE HP 59501B IS 706. ALL OF THE 
PROCEDURES MUST AGREE WITH THE FORTRAN77 PROGRAM THAT CALLS 


THEM WITH RESPECT TO THE NUMBER AND ORDER OF THE VARIABLES 
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ANO/OR STRINGS PASSED FROM ONE TO THE OTHER.) 


IMPLEMENT 


{ 


FOR MORE INFORMATION ON ALL THE I/O COMMANDS DEFINED BELOW, 
REFER TO THE PASCAL 3.0 PROCEDURE LIBRARY MANUAL; ESPECIALLY 
THE FOLD-OUT PAGES OF THE REFERENCE SECTION (P. 305).} 
PROCEDURE MY_WRITESTRINGLN3 
BEGIN 
WRITESTRINGLN(DEV;STR) $ 
WRITESTRINGLN OUTPUTS THE STRING; A CARRIAGE RETURN; AND A 
LINE FEED.} 
END; 
PROCEDURE MY_TRIGGERV; 
THIS PROCEDURE IS NOT USED IN THE FINAL VERSIONS OF THE THERMISTOR 
THERMAL CONDUCTIVITY TEST PROGRAMS BECAUSE THE THERMISTOR'S 
VOLTAGE IS READ AT THE EXTERNAL PORT OF THE DVM. THIS PROCEDURE 
WAS UTILIZED WHEN THE THERMISTOR WAS CONNECTED TO CHANNEL 09.} 
VAR 
FMT: STRING(80]; 
BEGIN 
FMT: ='AIO9';$ 
AI IS A COMMAND RECOGNIZED BY THE HP 34¢97A. AI, WHICH STANDS FOR 


ANALOG INPUT; CLOSES THE DESIRED CHANNEL AND TRIGGERS THE DVM TO 


TAKE A MEASUREMENT.) 

WRITESTRINGLN(DEV,FMT) 3 

END; 

PROCEDURE MY_WRITESTRING} 
THE MOST IMPORTANT PROCEDURE. THIS PROCEDURE BUILDS THE DIGITAL 
WORD THAT REPRESENT THE DESIRED CURRENT TO BE OUTPUTTED BY THE 
COUPLED HP 59501B AND HP 6186C. THE WORD BUILDING CONSISTS OF 
TAKING A TYPICAL NUMBER LIKE 2999 AND BUILDING A STRING BY 
ASSIGNING, IN ORDER AND ONE AT A TIME, THE 2, THE 9, THE 9) AND 
THE LAST 9.) | 

VAR 


FMT: STRING(80]; 
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{ 


ST: INTEGER; 
NX: INTEGERS 


RESULT: INTEGER; 


BEGIN ON 
FMT: = ° Uh 
ST: = 13 


ST DESIGNATES THE FIRST POSITION IN THE STRING FMT.} 

RESULT: = A DIV 1000; 
BY USING DIV, RESULT IS EQUAL TO THE FIRST NUMBER OF THE DIGITAL 
WORD (EG., 2)} 

STRWRITECFMT,ST,»NX,;RESULT: 1)3 
STRWRITE IS A TYPICAL PASCAL COMMAND THAT WRITES THE 2 INTO THE 
DESIGNATED POSITION WITHIN FMT. REFER TO PAGE 202 OF THE HP 
PASCAL LANGUAGE REFERENCE FOR MORE INFORMATION. } 

ST: = NX3 

A: = A-RESULT¥100035 
THE PREVIOUS CALCULATION GETS RID OF THE LEADING NUMBER (EG.; 
2999 BECOMES 999) AND THE PROCESS STARTS AGAIN.} 

RESULT: = A DIV 1003 

STRWRITECFMT»ST»NXsRESULT:1)3 

ST: = NX3 

A: = A-RESULT#®1003; 

RESULT: = A DIV 10; 

STRWRITECFMT,ST,»NX,;RESULT: 1)35 

ST: = NX3 

A: = A-~RESULT#10; 

RESULT: = A DIV 13 

STRWRITECFMT,ST;NX,RESULT: 1)3 

WRITESTRING(DEV, FMT); 
ONCE THE STRING FMT IS BUILT, WRITESTRING SENDS THE SPECIFIED 
STRING TO THE SPECIFIED DEVICE.} 

END; 

PROCEDURE MY_READNUMBER; 


BEGIN 
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READNUMBER (DEV, A) 
{ READNUMBER PERFORMS A FREE FIELD NUMERIC ENTRY FROM THE SPECIFIED 
DEVICE.) 
END; 
PROCEDURE MY_CLEAR; 
BEGIN 
CLEAR(DEV) ; 
{ CLEAR ATTEMPTS TO SEND A FORM OF THE CLEAR MESSAGE TO THE 
SPECIFIED DEVICE(S).) 
END; 


END. 
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200 


220 


230 


240 


250 


260 


270 


280 


290 


300 


310 


320 


330 


340 


THIS PROGRAM TAKES THE ASCII DATA FILE DEVELOPED FOLLOWING 
THE V¥¥2 VERSUS 1/(t#*.5) CALIBRATION STEP AND CONVERTS IT 
INTO A BOAT FILE THAT IS READILY UPLOADED INTO THE BASIC 


STATISTICS PACKAGE. 


THIS PROGRAM IS STORED UNDER THE FILENAME: CONVERSION. 


THIS PROGRAM WAS WRITTEN USING HP BASIC 3.0. 


DIM X(2200),XxX(2200) 

PRINT “ AT THIS TIME, PLACE THE DATA DISK WITH THE ASCII SOURCE DATA FILE” 

PRINT " IN THE RIGHT-HAND DISC DRIVE, #3." 

PRINT =>" 

PRINT “ ALSO, PLACE THE THE DATA DISC WHICH WILL STORE THE BDAT EQUIVALENT” 

PRINT " DATA FILE IN THE LEFT-HAND DRIVE, #4." 

PRINTS! 

PRINT | = 

PRINT " PRESS <CONTINUE> WHEN READY" 

PAUSE 

PRINT i = 

PRINTS es 

PRINT " " 

PRINT 5 

PRINT " WOULD YOU LIKE TO CATALOG EITHER DISK TO SEE WANT FILES PRESENTLY EXIST?” 
RAGE ENTER EITHER THE DESIRED DISC DRIVE NUMBER (3 OF 4) OR O TO CONTINUE. ",NOR 
PRINT See 


PRINTS Ss 


CATALOGING THE CONTENTS OF THE DISC IS DRIVE #3. 
IF NOR=3 THEN 

CAT ": INTERNAL" 

PRINT " " 

PRINT " " 


PRIN lene 
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350 
360 
370 
380 
390 
400 
410 
420 
430 
440 
G50 
460 
470 
480 
490 
500 
510 
520 
530 
540 
550 
560 
570 
580 
590 
600 
610 
620 
630 
640 
650 
660 
670 


680 


GOTO 240 


ENO IF 


CATALOGING THE CONTENTS OF THE DISC IN DRIVE #4. 


IF NOR=4 THEN 


CAT “: 
PRINT 
PRINT 
PRINT 
GOTO 2 
END IF 
PRINT 
PRINT 


PRINT 


INTERNAL 5431" 


G0 


PRINTING THE 12TH ASCII CHARACTER, A FORM FEED --> ACTS LIKE A 


CLEAR SCREEN COMMAND. 


PRINT CHR$(12) 


INPUT 
PRINT 
PRINT 
PRINT 
PRINT 


INPUT 


“ENTER THE FILENAME OF THE ASCII SOURCE DATA FILE (EX. 9CAS_1A__) 


“ENTER THE FILENAME OF THE BDAT EQUIVALENT DATA FILE" 
"(SUGGESTION: USE THE SAME FILENAME AS ABOVE, IE. "3M$3")." 
N$ 


DESIGNATING DRIVE #4 AS THE DEFAULT DRIVE. 


MASS S$ 


TORAGE IS “: INTERNAL;4,1” 


ALLOCATING SPACE ON THE DATA DISC IN ORIVE #4 FOR A NEW FILE. 


CREATE 


BDAT N$;70,256 


HAVE TO GIVE EACH I/O FILE NAME A CORRESPONDING PATH NAME. 


ASSIGN @PATH1 TO N$ 
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690 
700 
710 
720 
730 
740 
750 
760 
770 
780 
790 
800 
810 
820 
830 
840 
850 
860 
870 
880 
890 
900 
910 
920 
930 
940 
950 
960 
970 
980 
990 
1000 
1010 


1020 


DESIGNATING ORIVE #3 AS THE DEFAULT ORIVE. 
MASS STORAGE IS “: INTERNAL” 


ASSIGN @PATH2 TO M$ 


THE FIRST DATA VALUE (NEND) IS THE TOTAL NUMBER OF DATA VALUES 
IN THE ASCII FILE. 
ENTERING IN THE DATA FROM THE DATA FILE WITH THE CORRESPONDING 
PATH NAMED PATH2. 


ENTER OPATH2$NEND 


READING THE DATA FROM THE SOURCE (ASCII) DATA FILE 


PRINT CHRS$(12) 

PRINT ™ READING . . . . APPROXIMATELY A MINUTE" 
FOR I=1 TO NEND 

ENTER OPATH2;X(I) 


NEXT I 


CLOSING THE I/O PATH. 
ASSIGN @PATH2 TO * 
CHR$(12) IS THE ASCII CHARACTER FOR A FORM FEED -- CLEARS THE SCREEN 


PRINT CHRS$(12) 


WRITING THE DATA TO THE NEW (BDAT) FILE. 
PRINT “ WRITING .... ” 
MASS STORAGE IS “: INTERNAL;4;1" 

OUTPUTTING THE ARRAY OF DATA ALL AT ONCE. 


OUTPUT @PATH13 X(*) 


CLOSING THE I/0 PATH. 


ASSIGN @PATH1 TO * 


PRINT CHR$(12) 


PRINT ae 
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1040 PRINT " Lsholoiaial FINISHED Lslahaiaialt 

1050 PRINT => 

1060 PRINT “ " 

1070 

1080 PRINTING INFORMATION THAT FURTHER AIDS THE DATA REDUCTION SEQUENCE. 


1090 NN=NEND/3 


1100 FOR I= 
1110 PRINT 
1120 PRINT 
1130 PRINT 
1140 PRINT 
1150 PRINT 
1160 PRINT 
1170 PRINT 
1180 PRINT 
1190 PRINT 


1200 


1 TO 2 


USING 19A,2X;10A;"BDAT DATA FILENAME: “,N$ 


“THREE VARIABLES WERE STORED: V¥#2 (V¥#2), T_SELF-HEAT (C), TIME (SEC)" 


“TOTAL NUMBER OF OBSERVATIONS (FOR STATISTICAL DATA REDUCTION) = “3NN 


1210 PRINTING THE PERTINENT DATA TO THE HP PRINTER (701) AND THE CRT (1). 


1220 PRINTE 
1230 NEXT I 
1240 

1250 PRINTE 
1260 MASS S 


1270 END 
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